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Abstract 

Extending investigations of Metivier&Zumbrun in the hyperbolic 
case, we treat stability of viscous shock and boundary layers for viscous 
perturbations of multidimensional hyperbolic systems with characteris- 
tics of variable multiplicity, specifically the construction of symmetriz- 
ers in the low-frequency regime where variable multiplicity plays a 
role. At the same time, we extend the boundary-layer theory to "real" 
or partially parabolic viscosities, Neumann or mixed-type parabolic 
boundary conditions, and systems with nonconservative form, in ad- 
dition proving a more fundamental version of the Zumbrun-Scrrc- 
Roussct theorem, valid for variable multiplicities, characterizing the 
limiting hyperbolic system and boundary conditions as a nonsingular 
limit of a reduced viscous system. The new effects of viscosity are seen 
to be surprisingly subtle; in particular, viscous coupling of crossing 
hyperbolic modes may induce a destabilizing effect. We illustrate the 
theory with applications to magnetohydrodynamics. 
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1 Introduction 



This work is motivated by the stability analysis of boundary value problems 
and shock waves for viscous perturbations of multidimensional systems of 
conservation laws. In this analysis, three main steps are present. In the 
first step, one constructs simple waves w((u ■ x — at)/e), or "profiles" , which 
are exact solutions of the viscous equation, with viscosity of order e. This 
amounts to solving an ordinary differential equation (the profile equation); 
the solutions describe the fast transition between the hyperbolic solution and 
the parabolic boundary conditions (boundary layers) or between two smooth 
hyperbolic solutions (shock layers). Next, given a profile, formal plane wave 
or spectral analysis yields necessary stability conditions in terms of Evans 
functions. The second step is to compute explicitly this function on specific 
examples and check the stability conditions. The third step is to prove the 
linear and nonlinear stability of solutions, assuming that the suitable Evans- 
Lopatinski condition is satisfied, in particular for curved fronts or boundaries 
and non-piecewise constant hyperbolic solutions. This paper deals with the 
third step, with specific applications to magneto- hydrodynamics. The first 
and second steps are discussed for shock and boundary layers, respectively, 
in companion papers GMWZ6 and GMWZ5 . 

We concentrate on the construction of symmetrizers for the linearized 
equations, and more specifically in the so called low-frequency regime, as 
they are the key point in the proof of stability estimates which eventu- 
ally yield short-time existence and nonlinear stability theorems; see |Majl 
|Maj2| for hyperbolic shocks and Me ZiTTl IHM W X2I IHMWZ31 KtM W/4] for 
viscous perturbations. In these papers, it is proved that strong stability 
estimates hold, under the natural uniform Lopatinski condition, or Evans' 
condition, provided that the equations satisfy a structural condition, called 
the block structure condition (see Majl} IMaOsj in the hyperbolic case and 



[MeZnlj for the viscous case). This condition is in some sense necessary for 
the construction of Kreiss' symmetrizers which are used to prove the stabil- 
ity estimates. It is satisfied in the case of inviscid Euler's equations of gas 
dynamics ( |Majl| ), but does not hold in other interesting examples such as 
the equations of magneto- hydrodynamics (MHD). So there is a real need for 
an extension of the analysis beyond the class of systems satisfying the block 
structure condition. This is done in McZu2 for hyperbolic systems and the 
main goal of this paper is to extend the analysis to viscous systems, in view 
of applications to MHD. 

We carry out in passing several other useful generalizations of the basic 
boundary-layer analyisis of |MeZulllMe3j . extending the theory to "real" or 
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partially parabolic viscosities, Neumann or mixed- type parabolic boundary 
conditions, and systems with nonconservative form. In addition, we prove 
a more fundamental version of the Zumbrun-Serre-Rousset theorem, valid 
for variable multiplicities, characterizing the limiting hyperbolic system and 
boundary conditions as a nonsingular limit of a reduced viscous system as 
frequency goes to zero. Extensions to the shock case are given in GMWZ6 . 
Consider boundary value problems for hyperbolic systems 

d 

(1.1) d t u + AjdjU 

i=i 

on {xd > 0} with boundary conditions on {x^ = 0} which is assumed to 
be noncharacteristic. The plane wave analysis of such system leads to con- 
sider ordinary differential system in z > 0, depending on Fourier-Laplace 
frequencies £ = (r, 77, 7), with 7 > 0, 

d-l 

(1.2) d z -H ((), HoiO-^-A^^iT + ^ + Y^VjAj) 
Viscous perturbations of (|l.lj) are systems of the form 

d d 

(1.3) d t u + ^2 A i d i u ~ e d i ( B j,kdku) 

j=l j,k=l 

with natural structural conditions which are recalled below. The low-frequency 
plane wave analysis of such systems lead to consider perturbations of (|l,2j) : 

(1.4) d g -H(C,p)), H((,p):=H (()+pH 1 ((,p) 

depending smoothly on an additional parameter p > (see section 2 below 
or IMeZulllZlmTillMe!] ). 

In this paper, our main concern is the construction of symmetrizers 
p) for (|1.4|) . The precise conditions we impose on E are given in Section 
01 In particular, we focus on smooth symmetrizers, as they serve as symbols 
for pseudodifferential symmetrizers in the variable coefficient analysis. 

When p = 0, Eo(C) = E(£,0) is a symmetrizer for Ho(C)- Such sym- 
metrizers were constructed first for strictly hyperbolic systems (|1,1|) by 
Kreiss ([Kre]) (see also jChPij l. Strict hyperbolicity is used at only one 
place: it implies that Hq can be put in a normal form, which is called block 
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structure in |Majl IMaOsj . Therefore Kreiss' construction of symmetrizers 
extends immediately to systems which satisfy this block structure condition. 
In |MeZu2| . it is proved that this condition is satisfied if and only if the 
symbol A(£) := °^ P-^P i s smoothly diagonalizable for £ ^ 0, recov- 

ering known examples such as Euler's equation of gas dynamics or Maxwell's 
equations. The second important result in |MeZu2| it that the construction 
of symmetrizers is extended to a class of symmetric systems which are not 
smoothly diagonalizable for all £ 7^ 0: we demand that the "bad" multiple 
modes are totally incoming or totally outgoing (see the Definitions 14.11 and 
14.31 below), and this applies to inviscid MHD. 

In the small viscosity case, the construction of symmetrizers is performed 
in jMeZulj . with application to the analysis of shocks in |GMWZlllGMWZ2l 
GM WZ3| IGMWZ4] . assuming that the eigenvalues of A(£) have constant 
multiplicity for £ 7^ 0. As mentioned above, this assumption rules out the 
case of MHD. 

The main objective of this paper is to start the analysis of (|1.3[) or (|1.4|) 
when the constant multiplicity assumption is relaxed and in particular to 
investigate the construction of symmetrizers. It turns out that the influence 
of the viscosity is much more subtle than expected near multiple modes. In 
some cases, it may induce destabilizing effects. Let us list several new phe- 
nomena which can occur when there are multiple modes with nonconstant 
multiplicity. 

• Smooth diagonalization of Aq(£) implies a smooth block reduction for 
flo(C)- The perturbation pH\ in general couples the different blocks asso- 
ciated to a multiple eigenvalues (and this occurs for MHD). If the crossing 
eigenvalues do not have the same behavior with respect to the boundary 
(typically if they are not all incoming nor all outgoing), the spectral nega- 
tive space E~(£,p) is not continuous (in general) at p = 0. This happens 
for slow shock waves in MHD. This phenomena is excluded when the eigen- 
values have constant multiplicities; see |MeZu3| (in this case, since cross- 
ing eigenvalues are equal, they have the same behavior with respect to the 
boundary). Recall from MeZu2j that the continuity of E~ is a necessary 
condition for Kreiss' construction of smooth symmetrizers, more precisely 
for the existence of what is called below, smooth K-families of symmetrizers. 
In any case, the discontinuity of E~ is a major difficulty in the construction 
of symmetrizers. 

• As a consequence of the previous phenomenon, the Evans function 
can be discontinuous at p = 0. In the shock problem, the usual Evans 
function is in every case singular at p = (see |ZuSe| ). but the remark 
applies to the modified (or desingularized) Evans function introduced in 
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GMWZ3l lUiyLWZ4| (see below). 

• Because of the lack of continuity of the Evans function, it may happen 
that the strong Lopatinsky stability condition for the hyperbolic problem ( at 
p = 0) does not imply the strong Evans stability condition for small p. This is 
in sharp contrast with the known results obtained in the constant multiplicty 
case flZuSa IZumll iRoul IMeZull IGMWZ31 IOMWZ4Q . This is illustrated 
by an example in Section [3and this can occur for MHD, for some ad hoc 
boundary condition. An interesting question is to know whether this can 
happen or not for physical boundary conditions, in particular for slow MHD 
shocks. 

On the other hand, we prove in this paper the existence of smooth sym- 
metrizers under a natural generalized block structure condition for (|1.4|) . We 
also provide a geometrical characterization of this condition on the matrices 
A and B occurring in ()1.3j) . Moreover, modes that are totally incoming or 
totally outgoing do not cause trouble in the analysis of E~ nor of the Evans 
function. They are easily handled in the case of symmetric systems as in 
|MeZu2j . For instance, an important outcome of the present paper is the 
following result. We refer to the next sections for precise definitions. 

Theorem 1.1. Suppose that the full system ()1.3j) is symmetric. Suppose 
in addition that the eigenvalues of the hyperbolic system Ijl.lj) are either 
semi-simple with constant multiplicity or totally nonglancing in the sense of 
Definition \4.5\ Then, there are K-families of symmetrizers for the associated 
reduced system (|1.4|) . for p > sufficiently small. 

As recalled in the next section, K-families of symmetrizers provide Kreiss 
symmetrizers for boundary value problems which satisfy a uniform Lopatin- 
ski stability condition. One important application and motivation is the 
following 

Example 1.2. Fast Lax' shocks for MHD satisfy the assumptions of Theo- 
rem Al.li 

But we also have the following 

Counter example 1.3. Slow Lax' shocks for MHD do not satisfy the as- 
sumptions of Theorem \l.l\ 

Remark 1.4. Fast shocks with small magnetic field are perturbations of 
acoustic shocks of gas dynamics, whose stability has been studied by A.Majda 
( Majl ). Therefore, there are good reasons to think that the uniform Evans- 
Lopatinski condition is satisfied for Fast Lax' shocks for MHD, at least for 
perfect gases state laws and small magnetic field. 
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Remark 1.5. When the assumptions of Theorem 11,11 are not satisfied, or 
more generally when the generalized block structure fails, one could try to 
construct nonsmooth symmetrizers as in |MeZu2j . The counterexample ll.3l 
would be a good motivation for that. However, nonsmooth symmetrizer 
would require much more sophisticated pseudodifferential tools to handle 
variable coefficients. Moreover, slow shocks are not so closely related to 
acoustic shocks, and it is not known whether the uniform Lopatinski condi- 
tion is likely to be satisfied or not. 

When all the eigenvalues have constant multiplicity, Theorem 11.11 is 
proved in |MeZulj (see also |Me3| ) . 1 The construction is based on a re- 
duction of 1)1.4(1 to a suitable block diagonal form. Blocks which correspond 
to totally nonglancing modes (incoming or outgoing) are treated using the 
symmetry of the system as in Me Zu2| . For other blocks, we discuss in detail 
in Section 4 the generalized block structure condition which is needed for the 
construction of Kreiss symmetrizers. 

The symmetrizers are used in jMeZull IGMWZ11 IGMWZ21 IGMWZ31 
GMWZ4 to prove maximal stability estimates for boundary value problems. 
The Fourier multipliers S(p, £) serve as symbols for pseudo-differential sym- 
metrizers. All the other steps in these papers, linearization, paralineariza- 
tion, separation of frequencies, the high- and medium-frequency analysis, the 
conversion of the plane wave or symbolic calculus into an operator calculus 
via the use of a paradifferential calculus, are independent of the constant 
multiplicity assumption which was assumed there as a sufficient condition 
for the generalized block structure condition. Therefore, all these analyses 
are valid under the assumptions of Theorem 13.71 

As already mentioned, the main novelty of this paper with respect to 
previous works of the authors is the consideration of systems with variable 
multiplicity. To lighten the presentation, we will now we focus on bound- 
ary layers for noncharacteristic boundary value problems. The extension 
to classical, conservative Lax-type shocks requires only to incorporate the 
ideas already explained in detail in (for instance) | GMWZ3 , GMWZ4 . (The 
treatment of nonconservative and or undercompressive shocks involves new 
issues, and is carried out in [(3MWZ6 .) Similarly, we will concentrate only 
on the symbolic analysis for constant- coeffient equations and the construc- 
tion of smooth Fourier-Laplace multipliers. The passage from these multi- 
pliers to linear and nonlinear stability estimates for variable coefficients is 

lr The reduction to 11. 41 is carried out for strictly parabolic viscosities in MeZul Mc3 
and for partial viscosities in K1T1W 7.1} . However, the form of Hi is the same in each case 
(a consequence of Kawashima's genuine coupling condition, Assumption (H4) below). 
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already performed in previous works (see [MeZull IGMWZ31 IGMWZ4| ) and 

can be used as an independent black box. 

2 Spectral stability 

In this section, we recall the main steps of the spectral stability analysis of 
noncharacteristic boundary layers, refereeing to |Gue| IGrGu| IMeZul| IMe3| 
IZumll IZu2l IGMWZ3 , GMWZ4 for details and further references and appli- 
cations to the similar analysis of shock profiles. In particular, we give a new 
proof of the Zumbrun-Serre Lemma ( Zu Sel fRouj ) which allows for variable 
multiplicities. Moreover, not only does it provide a comparison between the 
Evans function of the viscous equation and the Lopatinski determinant of 
the inviscid system, but it also shows the link between the equations them- 
selves: for low frequencies, the viscous boundary value problem decouples 
into two boundary value problems, one of them being a nonsingular per- 
turbation of the limiting hyperbolic boundary value problem. We will also 
recall from [GMWZ4] the main arguments for the high-frequency regime. 

2.1 Structural assumptions 

Consider a system of equations 

d d 
(2.1) C e (u) := A (u)d t u + y £ j A j (u)d j u -ej = 0. 

j=l j,k=l 

When e = 0, £o is first order and assumed to hyperbolic; e plays the role 
of a non-dimensional viscosity and for e > 0, the system is assumed to 
be parabolic or at least partially parabolic. Classical examples are the 
Navier-Stokes equations of gas dynamics, or the equations of magneto- 
hydrodynamics (MHD). 

The form of the equations is preserved under a change of unknowns 
u = $>(u) or multiplication on the left by a constant invertible matrix. To 
cover the case of partial viscosity and motivated by the examples of Navier- 
Stokes equations and MHD, we make the following assumption: 

Assumption 2.1. (HO) The matrices Aj and Bj^ are C°° N x N real 
matrices of the variable u £ IA* C M. N . Moreover, for all the matrix 

Aq(u) is invertible. 

(HI) Possibly after a change of variables u and multiplying the system on 
the left by an invertible constant- coefficient matrix, there is N' G {!,..., N} 



S 



and there are coordinates u = (u l ,u 2 ) G R N N x M> N and f = (J 1 ,/ 2 ) G 
~ N x such that the following block structure condition is satisfied : 



(2.2) 4>(u) :=#(«) = jL), B iifc (u) = ( 





Bf k )> 



We refer to [GMWZ4j or |Zu2j for further comments and explanations. 
From now on we work with variables u = (u , u 2 ) G U* such that 1)2.2(1 holds. 
We set 

(2.3) Aj = fj, Aj = Aq Aj, Bj^ = Aq Bj^, 

and systematically use the notation M 01 ^ for the sub-blocks of a matrix M 
corresponding to the splitting u = (u ,u 2 ). Note that 



(2.4) B jM {u) := A Q (u)- L B jk (u 





n 



o bS(»), 



The triangular form of the equations also reveals the importance of the (1,1) 
block which plays a special role in the analysis : 



d 

(2.5) L ll {u,d) =Y J A f(u)d j , or L n (u,d) = {A 1 J{u))~ 1 L 11 (u,d). 

3=0 

In this spirit, the high-frequency principal part of the equation is 
(2.6) 



L U (u, d)u x 
d t u 2 - eB 22 (u,d)u 2 



with B 22 (u,^) = Yl'j k=i£j£kB 2 'Jl(u). We refer to Lemma [7.31 for a more 
detailed account of this notion of principal part. The first natural hypothesis 
is that L 11 (it, d) is hyperbolic and dt — B (u, d) is parabolic in the direction 
dt. 

Assumption 2.2. (H2) There is c > such that for all u €U* and £ G R d , 

the eigenvalues of B (u, £) satisfy Re/i > c|£| 2 . 

(H3) For all u£U* and all £ G R d \{0}, A U (u,£) = E^i&ApO) ^ 
onh/ real eigenvalues. 

For the applications we have in mind such as Navier-Stokes and MHD, 
the operator L is a transport field and (H3) is trivially satisfied. 
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Next we assume that the inviscid equations are hyperbolic and that 
Kawashima' s genuine coupling condition is satisfied for u, in some open 
sub domain U CU*. Let 

d d 

(2.7) A(u,0 = Y,^M U ) and B(u,0 = &kB j>k (u). 

3=1 j,k=l 

Assumption 2.3. (H4) There is c > such that for u £ U and ^ 6 W 1 , the 

eigenvalues o/M(w, £) + B(u,£) satisfy 

l£l 2 

(2.8) Ren > c- 1 



1 + 



Remark 2.4. (H4) implies hyperbolicity of the inviscid equation: for all 
u £ U and £ E M rf \{0} the eigenvalues of A(u,£,) are real. The set U may 
be thought of as the "hyperbolic set" where interior, inviscid solutions are 
be constructed, and the larger U* as the "hyperbolic-parabolic" set where 
exterior, boundary layer solutions are to be constructed, matching IA to 
boundary values in a multi-scale expansion. In contrast with GMWZ4] and 
|Zu2j , we do not assume here that the eigenvalues of A have constant multi- 
plicity. It is precisely the aim of this paper to substitute weaker conditions, 
allowing us to treat the case of MHD. 

Symmetric systems play an important role, and symmetry will be an 
important assumption in some of our results. In particular, the Assump- 
tion (H4) is satisfied when the following conditions are satisfied (see |Kawl 
EiS2j): 

Definition 2.5. The system (|2,1|) is said to be symmetric dissipative if there 
exists a real matrix S(u), which depends smoothly onu^lA, such that for all 
u G U and all £ G M rf \{0} ; the matrix S(u)Aq(u) is symmetric definite pos- 
itive, S(u)A(u,£) is symmetric and the symmetric matrix Re S(u)B(u,xi) 
is non negative with kernel of dimension N — N' . 

We consider a boundary value problem for (|2,1|) and the model case of a 
half space, which is given by {x > 0}, in some coordinates (j/i,... ,Vd-i, x ) 
for the space variables. We assume that the boundary is not characteristic 
both for the viscous and the inviscid equations. The principal term of the 
viscous equation is block diagonal as indicated in (|2.6|) The B 22 block is 
noncharacteristic by (H2). Restricting IA* to a component where the profiles 
will take their values, the condition for the A block reads 
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Assumption 2.6. U* is connected and for all u 6 U* , detA^ 1 (n) 7^ 0. 

For the inviscid equation, restricting U to the component where the 
hyperbolic solutions will take their value, the condition reads 

Assumption 2.7. U is connected and for all weW, det (Ad(u)) 7^ 0. 

By Assumption (H3) and Remark 12,41 A d (u) and Ad(u) have only real 
eigenvalues, which by Assumptions l2~71 and l2~o1 never vanish. This leads to 
two important indices : 

Notations 2.8. With assumptions as above, N+ denotes the number of 
positive eigenvalues of A^{u) for u G U and N\ the number of positive 

eigenvalues of (u) for u <EU*. We also set N b = N' + N]_ . 

The block structure ()2.6|) suggests that Nf, is the correct number of 
boundary conditions for the well posedness of (|2.1|) . for solutions with values 
in U* . Indeed, the high-frequency decoupling (|2.6|) suggests N' boundary 
conditions for u 2 and N+ boundary conditions for it . On the other hand, 
N + is the correct number of boundary conditions for the inviscid equation 
for solutions with values in U. Thus we supplement (|2.1|) with boundary 
conditions 

(2.9) T(u, edyU 2 ,ed x u 2 )\ x=Q = 0. 

Without pretending to maximal generality, we assume that they decouple 
into zero-order boundary conditions for u 1 and zero-order and first-order 
conditions for u 2 : 

( T 1 (n 1 ), ;E=0 = 0, 

(2.10) I T 2 (u% =0 = 0, 

[ T 3 (u,edyU 2 ,ed x u 2 )i x=0 = 0. 

with 

d-l 

T 3 (u,d y u 2 ,d x u 2 ) = K d d x u 2 + ^^Kj(u)djU 2 . 

i=i 

Assumption 2.9. Ti, T2 and T3 are smooth functions of their argu- 
ments with values in M. N + , ~R N ~ N and M> N respectively, where N" G 
{0, 1, . . • , N'}. Moreover, has maximal rank N" and for all uGW* the 
Jacobian matrices T'^u 1 ) and T' 2 (u 2 ) have maximal rank and N' — N" 
respectively. 
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2.2 Profiles and inviscid boundary conditions 

To match constant solutions u of the inviscid problem to solutions satisfying 
the boundary conditions, one looks for exact solutions of (|2.1j) (|2,9j) of the 
form: 

(2.11) u £ (t,y,x)=w(^j, 
such that 

(2.12) lim w(z) =u. 

The equation for w reads 

J A d {w)d z w - d z (B d4 (w)d z w) =0, z > 0, 
1 T(w,0,d z w 2 ) lz=0 = 0. 



(2.13) 



Solutions are called layer profiles. This equation can be written as a first 
order system for U = (w,d z w 2 ), which is nonsingular if and only if A\ is 
invertible (this indicates the strong link between Assumption 12.61 and the 
ansatz (f2~TTj) 1: 

(2.14) d z w 2 = w 3 , 

d z {B d>dW *) = {Af-Af{Af)^Af)w\ 

and the matrices are evaluated at w = (w ,w ). 

The natural limiting boundary conditions for the inviscid problem read 

(2.15) u\ x=0 e C, 

where C denotes the set of end points u such that there is a layer profile 
w G C°°(! + ;Zr) satisfying (l2~T2l (t2~T3l . The properties of the set C as well 
as the stability analysis of (|2.13|) depend on the spectral properties of the 
linearized equations from (|2.13|) near w(z). In particular we will discuss the 
notion of transversality for the profile w (see |MeZull IMe3j ). However, to 
avoid repetitions and prepare the multidimensional stability analysis, we en- 
large the framework and consider the multidimensional linearized equations 



from the full system (|2.1|) near solutions (|2.11|) . 

For further use, it is convenient to enlarge the class of functions w: 
consider a function C°°(M + ;U*) which converges at an exponential rate to 
and end state u £U: there is 5 > such that for all k G N 



(2.16) e Sz \d k z (w(z) -u)\ £ L°°(M + ). 
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We refer to such a function as a profile; it need not be a solution of (|2.13j) . 
though it will be in applications. Note that solutions of (|2.13[) 1)2.12(1 satisfy 
the exponential convergence above. 

Consider the linearized equations from (|2.1|) (|2.9j) around u £ = w(x/e): 

(2.17) C' u u = f, T'(u, sdyiL, ed x u)\ x = Q = g. 

Here T' is the differential of T at {w(0), 0, d z w(0)). C' Ue is a differential 
operator with coefficients that are smooth functions of z := x/e. Factoring 
out e~ l it also appears as an operator in edt,ed y ,ed x : 

(2.18) £' Ue = h(^,ed t ,ed y ,sd x y 

It has constant coefficients in (t, y), and following the usual theory of constant- 
coefficient evolution equations, one performs a Laplace- Fourier transform in 
(t, y), with frequency variables denoted by 7 + if and fj respectively, yielding 
the systems 

-L^,e(7 + if),iefj,ed x J. 

Next, we introduce explicitly the fast variable z = x/e, rescale the frequency 
variables as £ = (r, rj, 7) = e(f, fj, 7), and multiply the equation by e, reveal- 
ing the equation 

(2.19) L(z,-y + iT,ir),d z )u = f, T'(u 5 irju, d z u)\ z=0 = g, 



(2.20) L = -B(z)d 2 z + A(z, ()d z +M(z, (), 

with in particular, B{z) = Bd ) d{w{z)) and A u (z, (") = A^(w(z)). We do not 
give here the explicit form of A and M. Using (H2) and Assumption 12.21 
the equation is written as a first order system 

(2.21) d z U = G(z, C)U + F, T(C)U lz=0 = g, 
where 

(2.22) U = \u,d z u 2 ) = (u\u 2 ,d z u 2 ) G C N+N ' , 



(2.23) F=((A u (z))- 1 f\0,(B 22 (z))- 1 (-f 2 + A 21 (z)(A n (z)r 1 f 1 )). 
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The analysis of this equation depends on the size of the frequencies 
(. When £ is large, the character of the equations is dominated by the 
high-frequency principal part (|2.6[) . and we use a slowly- varying-coefficients 
analysis (related to the "tracking lemmas" of [Zumll [Zu2]) based on the 
relatively slow rate of change of coefficients compared to the size of the fre- 
quency; see |GMWZ31 IGMWZ4| and Section [7| below. For small or bounded 
frequencies £, we use the conjugation lemma of MeZul . The condition 
(|2.16|) implies that there is 5 > and an end state matrix G(u,C), depend- 
ing on the endstate u of w, such that 

(2.24) d k z (G(z,C)-G(u,()) = 0(e- 5z ). 

Lemma 2.10. Given £ G M a!+1 ; there is a smooth invertible matrix &(z, 
for z G M_|_ and Q in a neighborhood of (,, such that (j2,19|) is equivalent to 

(2.25) 8 Z U = G(u,()U + F, f(()U lz=0 = g. 

with U = $>(z,()U, F = $(z,C)F and f(C) = r(C)$(0,C). In addition, $ 
and converge the identity matrix at an exponential rate when z — > oo. 

Moreover, if the coefficients of the operator and w depend smoothly on 
extra parameters p (such as the end state u), then $ can also be chosen to 
depend smoothly on p, on a neighborhood of a given p. 

Remark 2.11. The linearized profile equations from (|2.13|) around w, are 
exactly (|2.19[) at the frequency £ = 0. In particular, Lemma I2.1UI implies 
that these equations are conjugated to constant-coefficient equations, via 
the conjugation by $(-,0). 

Next we investigate the spectral properties of the matrix G. Below, 
M.^ 1 denotes the open half space {( = (r, 77, 7) : 7 > 0} and M + + its closure 
{7 > 0}. We also introduce the matrices 

(2.26) P (u) := {B^)-\Af-Af{Af)^Af), 

d-l 

(2.27) H (u,0 := -(^(n))- 1 ((ir + 7)^oW + ^^A i (n)). 

j=i 

Lemma 2.12. i) Foru G li, Po(u) has no eigenvalue on the imaginary axis. 
We denote by the number of its eigenvalues in {Re/i < 0}. 

ii) For u S U and £ G M^ +1 \{0}, G(u,() has no eigenvalue on the 
imaginary axis. The number of its eigenvalues, counted with their multiplic- 
ity, in {Re fi < 0} is equal to N + + N 2 . = N b := N' + N]_ . 
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Hi) For a given u£U, there are smooth matrices V(u, Q on a neigh- 
borhood of (u, 0) such that 



(2.28) 



V- l GV 



H 
P 



with H(u, C) of dimension N x N, P(u, () of dimension N' x N', and 

a) the eigenvalues of P satisfy |Re/j| > c for some c > 0, 

b) there holds 

(2.29) H(u,C) = H (u,C) + O(\C\ 2 ) 

c ) at C = 0, V has a triangular form 



(2.30) 



V(u,0) 



Id V 
Id 



Proof, i) Take u G U. \iv 2 G kerP (^), then * (-(Af^Afv 2 , v 2 ) £ kerA d , 
implying that is not an eigenvalue of Pq. Similarly, if i£ is an eigenvalue 
of P then is an eigenvalue of i^A d + £ 2 -B^, which is impossible by (H4) if 
£ / is real. 

ii) Direct computations show that G(u,() = Gd{u, C) _1 M(-u, Q with 



with, in the splitting it = (u 1 ,^ 2 ). 
Sd(«) 



M 



®N'xN Mjv' x AT' 



On-N'xN' 

v E5i(u) , 



J = (On'xN-N 1 Un'xN') 



and 



i(u, C) = - J] iVj(Bj, d (u) + B dJ (u)) 



j'=i 

M(n, C) = (ir + 7)A)(«) + X] + X] VjVkBj^u) ■ 

3=1 j,k=l 

In particular, i£ is an eigenvalue of G(u, Q if and only if 7+ir is an eigenvalue 
of iA(rj,^) + B(r),£,), which, by (H4), implies either that 7 < if £ is real 
and (77, £) / or that £ = 0. 
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Thus G(u, C) has no eigenvalues on the imaginary axis and the number TV 
of eigenvalues in {Re/x < 0} is constant for u G IA and £ € M + +1 \{0}. That 
this number is equal to iV& = + is a consequence of the high-frequency 
analysis in Lemma 17.31 below (see also Lemma 1.7 in |Zu2| ). 

hi) Because M(u,0) = and A(u, 0) = Ai(it), there holds 



(2.31) G(u,Q) 



Since Po is invertible, G can be smoothly conjugated to a block diagonal 
matrix as in ()2.28j) . with V satisfying (|2.30f) and H(u,0) = 0. More precisely, 
the matrix V is 




(2.32) V 



( A d) A d P o l 



p -i 



The expansion (|2.29|) can be easily obtained by standard perturbation ex- 
pansions, and we refer to Lemma 14.231 below for a more precise version. 

For C, small, the number of eigenvalues of P in {Re/x < 0} is equal to 
N'i, and for 7 > 0, the number of eigenvalues of Hq(u,C) m the negative 
half space is constant, by hyperbolicity, and equal to N + . This implies that 
N = N + + Nl. □ 

Similarly, one considers the linearized equations from the inviscid hyper- 
bolic problem Cq{u) = around the constant solution u: 

(2.33) C'^ii = f. 

After performing a Laplace-Fourier transform, this equation reads 

(2.34) L (u,~f + iT,irj,d x )u = f 
or, with Hq defined at H2.21 



(2.35) 8 x u = H (u, C)« + A d \u)f. 

An important property for profiles is the notion of transversality (see 
MeZlTT] or [Me3] for the case of total viscosity). It concerns the linearized 
equations from (|2.11f) around w. As mentioned in Remark 12. 1H they cor- 
respond exactly to the first order system (|2,19|) with Q = 0. We abbreviate 
the homogeneous problem as 



(2.36) 



L(z,0,d z )w = 0, z>0, 
T(w,0,d z w 2 ) lz=0 = 0. 
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A corollary of Lemmas l2.1Ul and l2,12l is that the solutions of the homogeneous 
equation L(z,0,d z )w = form a space of dimension iV + N', parametrized 
by (u H ,u P ) G C N x C^' : 

(2.37) w(z) = <&h{z)u h + <S>p{z)e zPo ^u P 

where the matrices &h(z) and &p(z) are smooth and bounded on M + and 
$>h(z) — > Id as z — > Id. The solution is bounded if and only if up belongs to 
the negative space K~(Pq(u)) of Po(u), that is the invariant space of Pq(u) 
associated to the spectrum lying in {Re \i < 0}; thus the space S of bounded 
solutions has dimension N + N 2 . The space of solutions that tend to zero at 
infinity, denoted by So, has dimension iV£, corresponding to the conditions 
uh = and up G E~(Pq(u)). 

The boundary conditions in (|2.36|) read 

(2.38) T H u H + Tpup := T(w, d z w\ =Q = 0. 

Definition 2.13. The profile w is said to be transversal if 

i) there is no nontrivial solution w G So which satisfies the boundary 
conditions T(w, d z w 2 )\ z= o = 0, 

ii) the mapping w i— > T(w, d z w 2 )\ z= o from S to C Nb has rank Nf,. 

Equivalently, it means that ker T P nE~(Po(l*)) = {0} and that the rank 
of the matrix (T H ,T P ) from x E-(P (u)) to C^" is N b . 
If the profile satisfies condition i), there is a decomposition 

(2.39) = ¥ H ® F P , F P := TpE-(P {u)) 

with dimFf/ = iV + and dimFo,p = N 2 . Denote by tth and Tip the projec- 
tions associated to this splitting. 

For w G S given by (|2.37jl . one can eliminate up from the boundary 
conditions (|2.38|) and write them 

(2.40) Lred u H = 0, Up = R ,PU H , 

with 

(2.41) T red := ir H T H , R , P := -{T p y\ P T H 

and (£p) _1 is the inverse of the mapping T P from K~(Po(u)) to Fo,p. 

With these notations, ii) means that T_ Te d has rank iV_|_. Its kernel ker T red 
is the space of u G M rf such that there is a solution of w of (|2.36|) with end 
point u. It has dimension N — N + . 



17 



Remark 2.14. When uu is a layer profile, solution of (|2,13|) . the transver- 
sality condition implies that near the end point u, the set C in ()2.15|> which 
describes the limiting hyperbolic conditions is a smooth manifold of dimen- 
sion N_ = N — N + and kerr rerf is the tangent space to C at u. Therefore, 
the natural boundary condition for the linearized hyperbolic equation, and 
in particular for (|2,33|) . are 

(2.42) T red u = h. 

2.3 Evans functions and Lopatinski determinant 

For a given £ € M^ +1 \{0}, we now investigate the well-posedness of equation 
(|2~T3|) or equivalently (|2~2Tj) or (|2~23|) . Introduce the space E (£) of initial 
conditions (u(0), d z u 2 (0)) (or equivalently U (0)) such that the corresponding 
solution of L(z, £, d z )u = (or d z U — Q{z, QXJ = 0) is exponentially decaying 
at +oo. Lemmas 12.101 and 12 . 1 21 show that 

(2.43) E-(C) = *(0,C)E-(G(« > C)) 
where we use the following notations: 

Notations 2.15. Given a square matrix M, E _ (M) [resp. E + (M) denotes 
the invariant space of M associated to the spectrum of M contained in 
{Re n < 0} [resp {Refi > 0}]. 

In particular, by Lemma 12.121 E~(£) is a smooth vector bundle for £ £ 

I+ +1 \{0} and dim(E-(C)) = N b . 

The problems (|2.19() . 1)2.21(1 or (|2.25|) are well posed if and only if 

(2.44) E-(C)nkerT(C) = {0} or E~(G(u, ()) n ker f(C) = {0}. 

Note that, because the rank of T is at most iVj, and the dimension of E~ is 
Nb, this condition implies and is equivalent to 

(2.45) C N+N ' =E"(C)e ker T(() or C N+N ' = M~(G(u, ()) ker f(C). 
The Evans function is defined as 

(2.46) D(C) = \det N+N ,(E-(0MrT(0)\ 

where, for subspaces E and F of C n , det n (E, F) is equal to if dimE+dimF ^ 
n and is the n x n determinant formed by orthonormal bases in E and F if 
dimE + dimF = n. 
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Remark 2.16. The definition of the determinant above depends on choices 
of bases. Note that changing bases in E and F changes the determinant by 
a complex number of modulus one, thus leaves | det(E, F)| invariant. But it 
also depends on the choice of a scalar product on C n . Changing the scalar 
products (or changing of bases in C n )) changes the function det(E,F) to a 
new function det(E,F) such that c|det(E,F)| < |det(E,F)| < c^\ det(E, F)| 
where c > is independent of the spaces E and F. We will denote by 

(2.47) det » det or D^D 

this property. In particular, the definition of D is independent of the choice 
of orthonormal bases in E~ and ker T and all the uniform stability conditions 
stated below are independent of the choice of the scalar product. 

Remark 2.17. If the coefficients of the operator and the profile depend 
smoothly on parameters p, then the Evans function is also a smooth function 
of the parameters. 

These notations being settled, the weak stability condition, which is a 
necessary condition for well posedness in Sobolev spaces of (|2.17|) . reads: 

Definition 2.18. Given a profile w, the linearized equation (j2.17|) satisfies 
the weak spectral stability condition if D(Q) ^ for all Q G M + + \{0}- 

The next lemma is useful and elementary. 

Lemma 2.19. Suppose that E C C™ and T : C n i— ► C m , with rankT = 
dimE = m. If \ det(E, kerT)| > c > 0, then there is C , which depends only 
on c and |r*(rr*)~ 1 | such that 

VC/GE \u\<c\tu\. 

Conversely, if this estimate is satisfied then | det(E,ker T)| > c where c > 
depends only on C and \T\. 

Proof. Let it = r*(IT*) _1 r denote the orthogonal projector on (kerr) -1 -. 
Diagonalizing the hermitian form (ire, ire), yields orthonormal bases {ej} 
and {fj} in E and (kerT)- 1 - respectively, such that irej = Xjfj with < 
Xj < 1. Take any basis {g^} of ker T. Expressing the ej in the base {fk,gi}, 
implies that | det(E, ker T)| = Y\ Xj. Since Xj < 1 for all j, if this determinant 
is larger than or equal to c > 0, then minAj > c and for all e G E 

C \e\ < |vre| < |r*(rr*) _1 | |re|. 

Conversely, if the estimate is satisfied, then |e| < C\T\ \ire\ since Ye = Tire 
for all e G E. Therefore AjC|r| > 1 and the determinant is at least equal to 
(C\T\)- m . □ 
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There are analogous definitions for the linearized hyperbolic problem 
(|2.33f) with boundary conditions (|2.42|) . For 7 > 0, Hq(u,£) has no eigen- 
values on the imaginary axis, as a consequence of the hyperbolicity as- 
sumption (see Remark I2.4|) . The Lopatinski determinant is defined for 
C G R d + +1 := {7 > 0} by 

(2.48) D Lop (() = I det (E-(# (u,C),kexr red )|. 

By homogeneity of Hq, this determinant is homogeneous of degree zero in £ 
and one can restrict attention to ( G S d = = 1}. 

Definition 2.20. The linearized equation (j2.33|) ()2.42j) satisfies the weak 
spectral stability condition if -Dlop(C) 7^ f or a ^ C G 

2.4 Uniform spectral stability and maximal estimates 

The weak stability conditions and the reduction to constant coefficients of 
Lemma 101 guarantee the well posedness of (j2.19[l for fixed ( G M + + \{0} 
and in particular estimates of the form 

(2.49) |M| L2 + \\d z u 2 \\ L 2 + K0)| + \d z u 2 (0)\ < C(C)(||/|U2 + \g\). 

The next step in the study of (|2.17|) , is to perform an inverse Fourier-Laplace 
transform and thus requires suitable estimates for the solutions of (|2.19jl . 
with a precise description of the constants in the estimate above. 

By continuity in the weak stability condition implies that the estimate 
()2.49j) is satisfied with a uniform constant C when £ remains in a compact 

subset of \{0}. Thus the true question is to get a detailed behavior of 
the estimate when Q — > and when |£| — > 00. 

2.4.1 Low and medium frequencies 

Consider first the low-frequency case. Following |MeZulj . the uniform sta- 
bility condition reads: 

Definition 2.21. Given a profile w, the uniform spectral stability condition 
for low frequencies is satisfied when there are c > and po > such that 
D(Q > c for all ( G with < |C| < Po- 

By Assumption the rank of T(£) is always Nf,, and the norms of r(£) 
and (rr*) -1 are uniformly bounded for £ bounded. Thus, by Lemma 12.191 
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the low-frequency uniform stability condition holds if and only if there are 
C and po > such that 

(2.50) VCG< +1 , 0<\(\<p , WeE-(C): |*7| < C\T(()U\. 

Following MeZul], the expected maximal estimates for low and medium 
frequencies for the solutions of (|2.19|) read 

VIMU 2 (K + ) + l|d 2 u 2 || L 2 (R+) + W(0\ + |<^u 2 (0)| < 
(2-51) A 

where ip = (7 + |C| 2 )5 with C independent of ( G 1+ +1 \{0}, |C| < Po- Note 
that for fixed |£| > 0, this estimate implies (|2.49|) . 

The estimates (|2,51j) correspond to estimates for the solutions of the first 
order system (|2.21(l : 

(2.52) vWtfWv^ + \\U 2 \\ L 2 (R+) + \U(0\ < C{^\\F\\ L 2 (R+) + \g\) 

where U = (U 1 , U 2 ) £ C N x C^'. For the constant-coefficient system (j2~23|) 
the expected estimates read : 

(2.53) vW&Wv^) + \\U 2 \\lHr + ) + \U(0\ < C(-||F|| i2(M+) + \g\). 

Lemma 2.22. The estimates (j2,53|) imply 1)2.53(1 which imply ()2.51|) . 

Proof. (See MeZul ). Clearly, (|2.51f) is a particular case of (|2.52() applied to 
source terms F of the special form (|2.23() . Moreover, using the conjugation 
Lemma l2.10l there holds U = 0(1)U and U = 0(1)U and similar estimates 
for F and F. Moreover, 

U 1 = 0(1)U, U 2 = 0{e- dz )U l + 0(1)U 2 

with 6 > 0. We use the inequality 

\\e- Bz U l \\ L *< \U l {p)\ + \\d,U 1 \\». 

Moreover, the form of G(u, £) at £ = shows that 

dzU 1 = 0(\(\)U 1 + 0{l)U 2 + F 1 . 

Therefore, 

ll^ 2 ll^ < \\u 2 \\ L 2 + 1^(0)1 + Idllf/ 1 !!^ + hf 1 !!^. 

Since |C| < ip, this shows that (f233~j) implies l(235|> . □ 
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For ( in a compact subset of M < J +1 \{0}, all these estimates are true under 
the weak stability condition (see e.g. [MeZulj ). Note also (taking / = 
in 1)2.51)1 ) that the uniform stability condition ()2.50)l is necessary for the 
validity of the maximal estimate. The main subject of this paper is to prove 
that the uniform stability condition implies the maximal estimate 1)2.51)1 for 
low frequencies, under structural assumptions on the system weaker than in 
jMeZn 1 1 imvTWZ3l imVTWZ4j . allowing for instance to consider MHD. 

2.4.2 High frequencies 

For the high-frequency analysis, the maximal extimates that are proven in 
GMWZ4 concern homogeneous boundary conditions (g = 0) and read 

(1 +7)||u 1 || L 2 (M+) +A||u 2 || L 2 (]R+) + \\d z u 2 \\ L 2 {R+) 

(2.54) +(l + 7 )^|n 1 (0)|+A2|n 2 (0)| + A'^\d z u 2 (0)\ < 

C'(||/ 1 ||L 2 (R + )+A- 1 ||/ 2 || i2(R+) ). 

where A is the natural parabolic weight 

(2-55) A (C) = (1 + t 2 + 7 2 + H 4 ) 1/4 . 

The balance between the weights for u 1 and for u 2 is subtle: these com- 
ponents are decoupled in the high-frequency principal system 1)2.6)1 and the 
weights depend on their actual coupling through the nondiagonal terms and 
the boundary conditions. Here we see the importance of the form ()2.1U)1 
of the boundary conditions. Their linearized version, T'(u,ir]u 2 ,d z u 2 ) = g 
reads 

f IV(0) :=r i («; 1 (0)).« 1 (0)=^ 1 , 

(2.56) <^ T 2 n 2 (0) := V 2 (w 2 (0)) ■ n 2 (0) = g 2 , 

( r 3 (C)(u 2 (0),^ 2 (0)) :=K d d z u 2 (0)+K tg ( V )u 2 (0)=g 3 . 

with 

d-1 

(2.57) ^t g = E^ if iW°)) 

The complete maximal estimate with nonvanishing boundary source terms 
g, reads 

(1 + 7)||n 1 || i 2 (M+) +A||-u 2 || L 2 {M+) + ||^u 2 || L 2 (R+) 
+(l+7)V(0)|+AV(0)|+A-5|d z u 2 (0)| < 

c(II/ 1 IIl 2 (r+) + a ~ II/ 2 IIl 2 (r+)) 

+ C((l + 7 ) Vl + A Vl + A-i| 5 3 |) 
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with C independent of ( E M + large. Taking / = 0, this implies the 
following necessary condition : there are C and p\ > such that 

VCe< + \ ICI>Pi, VC/ = (« x ,u 2 ,u 3 ) gE-(C) : 

(2.59) (1+7)51^1 +A5|u 2 )| +A-5|n 3 | < 

c((i + T )^|r 1 « 1 | + Al|r 2 ^ 2 )| + A-^|r 3 (C)(^ 2 ,« 3 )l) 

This can be reformulated in terms of a rescaled Evans function (see [M eZulj : 
In <C N+N and C Nb introduce the mappings 

J c (u\ U 2 ,u 3 ) := ((l + 7)i« 1 J A5u 2 ,A-iu 3 ) 

(2.60) j x , 

Jd9\9 2 ,9 3 ) ■= ((1 +7) 5 <?\A^ 2 ,A"9 3 )- 
Note that J ? r(C)t/ = r sc (C) J C C/ with 

(2.61) r sc c/= (riu^^iQ^ + A-^tg^y). 

Thus (12331) reads 

(2.62) V£7 G J C E-(C) : |C/| < C| J ? r(C) J ? _1 C/| 
Introducing the rescaled Evans function 

(2.63) D SC (C) = |det(J c E~(C),J c kerr(C))|. 

we see that this stability condition is equivalent to the following definition: 

Definition 2.23. Given a profile w, the linearized equation (j2.17l) satisfies 
the uniform spectral stability condition for high frequencies when there are 
c > and p\ > such that D sc (() > c for all ( € with \(\ > p\. 

Note that for £ in bounded sets, Jq and J^ 1 are uniformly bounded and 
-D(C) ~ D SC (Q, thus the condition D SC (Q ^ is nothing but a reformulation 
of the weak stability condition. 

By Lemma 12.191 the high-frequency uniform stability is equivalent to 
(|2.59|) . In section E[ we will recall from [GMWZ4 that the uniform spectral 
stability implies the high-frequency maximal estimates (j2.58|) . under struc- 
tural assumptions on the system that are satisfied in many examples, includ- 
ing Navier-Stokes and MHD. 

Remark 2.24. The structural assumptions we refer to are connected with 
well-posedness of the initial-value problem for the viscous equations. For 
shock waves, they by themselves guarantee spectral stability and maximal 
estimates lGMWZ4| . For boundary- value problems, they reduce spectral 
stability to well-posedness of the frozen-coefficient boundary- value problem 
at the boundary; see |MeZul| IGM WZB for further discussion. 
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2.4.3 The inviscid case 

There are analogous definitions for the linearized hyperbolic problem (|2.33|) 
with boundary conditions 1)2.42(1 . Recall that the Lopatinski determinant is 
defined at (|2.48|) . Definition 12 . 2UI of weak stability is strengthened as follows. 

Definition 2.25. The linearized equation ()2.33|) (|2.42j) satisfies the uniform 
spectral stability condition holds when there are c > such that -D(C) > c 
for all (eSi:=S d n{ 7 >0}. 

This uniform stability condition is equivalent to a uniform estimate for 
all C G Sf: 

(2.64) Vu E E~(H (u, ()) : |u| < C\T red u\ 

The expected maximal estimates for solutions of ()2.33|) ()2.42|) are 

(2.65) lHu\\ L 2 + | u (0)| < C( 7 ^||/|| L 2 + \h\) 
with C independent of C G 

2.5 The Zumbrun-Serre-Rousset Theorem and the reduced 
low-frequency problem 

In this section, we extend the previous results of [ZuSej and |Rouj which 
link the low-frequency uniform stability of the viscous regularizations and 
the uniform stability of the limiting inviscid problem. First, we recall that 
the transversality of the profile is a necessary condition. 

Proposition 2.26. Given a profile w, if the low-frequency uniform spectral 
stability condition is satisfied, then w is transversal. 

Proof. Lemma 12.121 implies that for £ ^ small enough, U is a solution of 
(t2~2T|) if and only if t (u H ,u P ) = V-^QU satisfies 

(2.66) d z u H = H(u,()u H + fH, 

(2.67) d z u P = P{u,()up + f P , 

(2.68) r H (C)uir(0)+rp«)up(0) := T(()U(0)=g, 

where *(/#, fp) = V~ 1 (()F and r# [resp Fp] denotes the restriction of rV 
to C N x {0 } [resp. {0} x C^']. In particular, 

E - (G(n, 0) = V(Q (E-(H(u, 0) E~(P(u, ())) . 
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With (|2.5U|) . this shows that the low- frequency uniform stability condition 
holds if and only if there are C and po > such that for all ( G with 

o < ici < po 

(2g9) Vuh €!"(#(«, 0), V«pGE-(P(u,C)): 

|«p| + |«p| < c|rp(c)«# + r P (c)up|. 

In particular, 

(2.70) \/u P £ E _ (P(u,C)) : |«p| < C|rp(C)up|. 



By Lemma |2.121 E~(P(u, Q) is a smooth bundle for ( in a neighborhood 
of 0. Moreover, r(£) and Tp(£) are smooth around the origin. This im- 
plies that \up\ < C\Tp(0)up\ on E _ (P(u, 0)), implying that condition i) of 
Definition 12.131 is satisfied. 

Since dim(E _ (G(C)) = rankf(£) = iVj,, implies that for all h G 

C^" and all ( G ij 1 " 1 with < |C| < po, there is 17(C) = V"(C)(«h(C),«p(0) 
in E-(C) C V(C)(C N ©E-(P(C))) such that f (C)f7(C) = ^ and |P(C))| < 
c|/i|. By compactness and continuity, letting C tend to zero, implies that 
there is U = V(0)(u H ,u P ) in V{0)(C N © E~(P(0))) such that f (0)£7 = h, 
showing that condition ii) of Definition 12.131 is also satisfied. □ 

Suppose that the profile w is transversal. Then, by i) of Definition 12.131 
and Remark 12.111 Tp(Q is an isomorphism from E~(P(u, £)) to its image 
Fo,p when £ = 0; by continuity this extends to a neighborhood of the ori- 
gin and the decomposition (|2.39|) valid at £ = 0, extends smoothly on a 
neighborhood of the origin: 



(2.71) C Nb = ¥ H © Fp(C), Fp(C) := T P (QM-(P(u, ()). 

Denote by ith(() arid vrp((") the projections associated to this splitting and 
define the reduced boundary operator as 

(2.72) T red (C) :=MC)rp(C) 5 
as well as the reduced boundary value problem 

(2.73) d z u H - H{u, C)u H = f H , r red (C)np(0) = h. 
The reduced Evans function is 

(2.74) D red (0 = | det {E~{H{u, C)),ker T red (0) |- 
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Definition 2.27. The reduced uniform stability condition is satisfied if 
Dred{C) > c > for all £ G R + \{0} wii/i |£| small enough. 

This is equivalent to the condition 

(2.75) Vu G E-(F(« ; 0) : |«| < C|r red (C)u|, 
for C G l d+1 \{0} small. 

Theorem 2.28. Given a profile w, the linearized equation (|2.19|) satisfies 
the low-frequency uniform spectral stability condition if and only if 

i) w is transversal, 

ii) the reduced problem (|2.7H|) satisfies the reduced uniform stability 
condition. 

Proof. We have already shown that the low-frequency uniform stability re- 
quires that w is transversal. Moreover, using the splitting (|2.71j) . we see 
that the uniform stability conditions (|2.5Uf) or (|2.69|) are equivalent to 

(2.76) \u H \ + \u P \ < c(\T red u H \ + \T P up + irpT H UH\j 

for all uh G K~(H) and up G E~(P) (to lighten notations we have omitted 
the (" dependance). Since Tp is surjective from E~(P) onto Fp, for all 
uh G K~(H) there is up G E~(P) such that Tpup = —ttpThuh and (|2,76|) 
implies (|2?T5|) . 

Conversely, if the profile is transverse, the estimate ()2.7(H) is valid at 
C = and extend by continuity to £ in a neighborhood of 0. With (|2.75|) . 
this clearly implies (|2,76|) , Q 

It remains to link the reduced uniform stability condition to the uniform 
(Lopatinski) stability condition for the hyperbolic boundary value problem, 
that is for the problem (|2.33f) with boundary conditions ()2.42|) . Note that 
these boundary conditions are given by T red = T re d(0) (see Remark l2.14|) . 

Because H vanishes at £ = 0, it is natural to use polar coordinates: 

(2-77) ( = pC, p=\C\, (£S d . 

In these coordinates 

(2.78) H(u,() = pH(u,(,p), H(u,(,p)=H (u,0 + O(p). 

Changing z to z = pz, u(z) to u{z) and f{z) to pf(z) the reduced problem 
(|2.73f) is equivalent to 

(2.79) d- z u H - H(u, C, p)u H = f H , r red (()u H (0) = h, 
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which, for p = 0, is exactly the inviscid problem (|2.35f) 1)2.421) . We are thus 
led to a nonsingular perturbation problem. 

Clearly, for (e^:=S J n{i> 0}, there holds E-(H(u,()) = 
E-(H(u,(,p) and D red (C) = D((,p) with 

(2.80) D(C, p) = | det {E-(H(u, (, p), ker T red (p()) \ 

Remark 2.29. For 7 > 0, Hq(u,C) nas no eigenvalues on the imaginary 
axis, as a consequence of hyperbolicity (see Remark 12. 4|) . By perturbation, 
this property holds true for H(u, C, p) for p small enough (depending on 7 > 
0). This shows that the vector bundle M~(H(u, C,p) which was defined on 
S + x]0, po] has a smooth extension to G S + x [0, po]> a s well as -D. Comparing 
with the definition of the Lopatinski determinant (|2.48|) . we see that 

(2.81) D Lop (C) = D((,0), for7>0. 

The next theorem, combined with Theorem l2.281 extends Rousset's the- 
orem |Rouj (see also [ZuSej for shocks). 

Theorem 2.30. Given a transverse profile w, if the reduced uniform spec- 
tral stability condition is satisfied, then the linearized hyperbolic problem 
(j'2.33f) ()2.42|) satisfies reduced uniform stability condition. 

Conversely, if the linearized hyperbolic problem is uniformly stable and 
the vector bundle E~ (H (u, ( , p) has a continous extension to S + x [0, po], 
then the reduced uniform spectral stability condition is satisfied and the lin- 
earized problem (|2.17j) satisfies the uniform low-frequency stability condition. 

Proof. The uniform estimate (|2,75|) implies that 

M < C\T red (()u\ 

for u G E~ (H (u, C , p) , C 6 S+ and p > small. If 7 > 0, every term is 
continuous up to p = and the estimate above implies (|2,64|) . that is 

|u| < C\T red {0)u\ 

for u € E _ (^o(M)C)) C ^ Si_. This implies that the hyperbolic problem in 
uniformly stable. 

If E~(H(u, p) has a continous extension to S + x [0, po], the reduced 

Evans function is has a continuous extension to S+ x [0, po]. The hyperbolic 
uniform stability and ()2.81j) imply that 

D(C,p) >c>0 
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for £ G S+ and p = 0. By continuity, this extends first to C G 5 + and next 
to p G [0, pi] for some p\ > 0. □ 

Remark 2.31. It is proved in |MeZu3| that when the eigenvalues of the hy- 
perbolic symbol A(u,t;) have constant multiplicity, and more generally when 
there is a smooth K family of symmetrizers (see the definition below), then 
the vector bundle E~(i?(u, Q,p) has a continuous extension to p = 0. The 
main concern of this paper is to construct K-families for systems with vari- 
able multiplicity. This is possible under suitable assumptions, and therefore 
the two theorems above extend a result of F.Rousset |Rou| . However, we 
will also show that the bundle E does not always admit a continuous exten- 
sion, with the result that the hyperbolic problem can be uniformly stable 
while the viscous problem is strongly unstable in the low-frequency regime. 
This seems to be a new phenomenon. 

Assuming transversality of w, Theorem 12.281 implies that the uniform 
spectral stability for low frequency is equivalent to the spectral stability for 
the reduced problem. There is an analogue for maximal estimates. The 
maximal estimates for the reduced problem (|2.79|) read 

(2.82) & + p)*\\&h\\v> + IM0)| < C[tf + p)-*\\f H \\v> + \h\) 

with C independent of £ G and p G]0, po]. Note that for p = and 7 > 0, 
this is the maximal estimate for the inviscid problem. Scaling back to the 
original variables, this estimate is equivalent to 

(2.83) (7 + |Cl¥lM| L2 + Wh(0)\ < C(( 7 + \C\ 2 r l i\\fH\\Li + \h\) 
for the solutions of ()2.73j) . 

Theorem 2.32. Suppose that the profile w is transversal. Then the maximal 
estimates (|2,53|) are valid for low frequencies if and only if the maximal 
estimates (|2.82|) for the reduced problem hold true. 

Proof. By Lemma \2.12\ P(u. C) has no purely imaginary eigenvalues. Thus, 
using symmetrizers (see e.g. [MeZulj and Section |31 below), there holds 

(2-84) H\\ L 2 + |«+(0)| < \\f+\\v, 

(2-85) ||up|U 2 < ||/p||^ + |«p(0)|, 

where ± denotes the smooth projections on the spaces E ± (P(?x, £). 
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The splitting (|2,71|) implies that the boundary condition (|2.68f) reads 

ir H g = T red u H (0) + ir H T P u^(0), 

ir P g = T P u p (0) + ir P T H u H (0) + 7TpTpu+(0). 

Moreover Tp is invertible on E~(P), hence |rp«p(0)| |«p(0)| and 

|r re dWfl-(o)| < \ir H g\ + |«p(o)|, 

I«p(0)| < \TTpg\ + \u H {0)\ + \u+(0)\. 
Suppose that the estimate (|2.83|) is satisfied. Then, 

(pWuhWl* + l«ff(°)l < ^ _1 II/h||l2 + \7r H g\ + \up(o)\. 

With (|2~HH) . this implies that 

(P\\u h \\l2 + \\ u p\\l2+\ u h(0)\ + hp (0)| 

<^ 1 II/h||l 2 + II/pIIl 2 + bl + 14(0)1- 

Thus, with (j2.84j) . we obtain that 

^II^IIl 2 + II^pIIl 2 + \uh(o)\ + 1 up (0)1 < ^ _1 ||/h||l 2 + II/pIIl 2 + M- 

Because V(u,0) has the special form (|2..30j) . f7 = V(uh,up) = (U 1 ,U 2 ) 
satisfies 

C/ 1 = 0(l)u H + 0(l)u P , U 1 = 0{\(\)u H + 0(l)u P 
Therefore, the solutions of (|2.25|) satisfy 

tpWU'Wv + \\U 2 \\ L2 + \U(0)\ < tp-^FWv + \g\. 
that is the maximal estimate (|2.53|) . 

Conversely, assume that the maximal estimate (|2.53|) is satisfied. Sup- 
pose that uh is a solution of (|2.66|) . By transversality, Tp is surjective from 
E~(P, () to its image Fp(£) and there exists there is up(0) in E - (P, Q such 
that 

(2.86) Tpu P (0) = -7T P r H u H (0) G F P (C). 

Consider up = e zP up(0) which is well defined and rapidly decaying at 
infinity since itp(0) £ E~(P, Q. It is a solution of (|2.67f) with fp = 0. 
Then U = V(uh,u p ) is a solution of QTEty with F = V{f H ,0). Thus 
(up, up) = V~ 1 U and there holds 

\Wh\\l* < WWl^ Kff(0)| < |f/(0)|, ||F|| L2 < \\f H \\ L 2 

and, by (l2~HH|) . f [7(0) = Tpup(O) + r P tip(0) = r red up(0). Thus the esti- 
mate 1)2 .53(1 immediately implies 1)2.83)1 . □ 
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3 Low frequency analysis: the main results 

This section is mainly devoted to the study of the reduced equation (|2.79f) , 
which is a nonsingular perturbation of the inviscid problem (|2.33|) , Our 
goal is to perform an analysis without assuming constant multiplicity of 
eigenvalues, thus allowing for examples such as MHD. The inviscid case is 
considered in MeZu2], and we want to extend the results to small viscous 
perturbations. 

3.1 Symmetrizers 

Consider the constant-coefficient linear first order system (|2.25[) . For clarity, 
we drop the tildes and reserve the notation u,U... for the unknowns and call 
p 6W the parameter called u in this equation, which now reads 

(3.1) d z U = G(p,()U + F, r(p,QU(0) = g. 

To prove energy estimates for the solutions of this equation, the main step is 
to construct symmetrizers. They are self adjoint matrices £) such that 

(3.2) Re(S(p,C)G(p,C)) >0. 

The symmetrizer is adapted to the boundary conditions and provides max- 
imal estimates for the traces when 

(3.3) S(p,C) > on kerT(p,C). 

The construction of such symmetrizers is in two steps: first, one con- 
structs a family of symmetrizers S K , which is independent of the boundary 
conditions; second one uses the uniform Lopatinski or Evans condition, to 
prove that if k is large enough then the simmetrizer is adapted to the bound- 
ary condition. 

More precisely, one considers a splitting 

(3.4) c^+^E-faOeE+fac) 

where E~ (p, Q is the negative invariant space of G(p, £) as above while 
E + (p, C) can De chosen arbitrarily so that the splitting (|3.4|) holds . De- 
noting by II (p, C) the projectors associate to this splitting, the family of 
symmetrizers S K is searched so that 

(3.5) s K > m(K)(n+)*n+ - (rr)*ir 
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where m{n) — ► +00 as k — > +00. 

Since E~ = kerll + , the stability condition (|2.44|) which reads 

(3.6) kerr(p,C)nE-(p,C) ={0} 
is also equivalent to an estimate 

(3.7) \U~u\ 2 <C{\ru\ 2 + \U + u\ 2 ). 

Therefore, if the family S K satisfies (|3.5[) . then for k large enough, there 
holds 

(3.8) s K > dd - c"r*r, c > 

and therefore S K is adapted to the boundary condition T. 

If Re T, K G > <5 K Id, then multiplying the equation by T, K and integrating 
by parts yields the estimate 

(3-9) SkWWI* + c\U(0)\ 2 < U\F\\ 2 L2 + C'\g\ 2 . 

Ok 

This is the sketch of the general argument. To obtain usable estimates, 
uniform versions of (|3.5|) (|3.7|) are needed as well as more precise versions of 
(J3.2|) (see below). Note that in this approach, the construction os symmetriz- 
ers is completely independent of the boundary conditions, and in particular 
of the validity of the stability conditions. In this paper, we concentrate on 
the construction of families of symmetrizers which satisfy (jHHj) . They are 
called K- families in }MeZu2j . 

3.2 Main results 

The construction of symmetrizers for middle frequencies, is performed in 
[MeZulj . By Lemma I2.12| the matrix G(p,C) has no eigenvalues on the 
imaginary axis when (£l \{0}. Therefore, 

Lemma 3.1. For all ( £ R \{0}, there is a neighborhood of (p, £) in 
U x such that for (p, Q in this neighborhood there is a smooth splitting 

(3.10) C"' =E-(p,C)®E + (p,C). 

where ^^(p, C) denote the invariant space of G(p,() associated to the spec- 
trum in {±Re/U > 0}. Denoting by H^(p, C) the smooth spectral projectors 
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associate to this splitting, there is a smooth family Y, K (p,() of self adjoint 
matrices such that for all (p, £) in the given neighborhood and all k > 1: 

i) Re T, K G > 0, 

( ' > ' 1M a) ReS K > K (n + )*n+ - (ir)*ir. 

Corollary 3.2. // the weak spectral stability condition is satisfied, then for 
all C G M \{0}, there are a constant C and a neighborhood of (p, £) in 
IA x M. d+1 such that for (p, Q in this neighborhood the solutions of (|3.1() 
satisfy 

(3-12) \\U\\ L 2 + \U(0)\ <C(\\F\\ L2 + \g\). 

We now concentrate on low frequencies. By Lemma 12.121 the matrix 
G(p, is locally smoothly conjugated to a block diagonal matrix (|2.28|) with 
diagonal blocks with H (p, £) of dimension N x N and P(p, C) OI dimension 
N' x N'. The system (|3.1j) is therefore equivalent to the equations (|2.6(ifl 
(j'2.67f) coupled by the boundary conditions (|2.68j) . 

In the block diagonal reduction ()2.28|) . we construct symmetrizers 

(3.13) V = ( S « ^ 

such that the properties (|3.2|) ande (|3.5j) are satisfied for each block inde- 
pendently. 

The construction of symmetrizers for the elliptic block P is standard 
and identical to the construction for middle frequencies, since P(p, 0) has 
no eigenvalues on the imaginary axis. Denote by Ep(j>, £) the subspaces 
of , invariant for P(p,(), associated to the spectrum in {±Re^ > 0}. 
Thus, for (p, £) in a neighborhood of (p, 0), there is a smooth splitting 

(3.14) C Ar '=E p 0E+ 

Denote by Tlp(p, C) the smooth spectral projectors associate to this splitting. 

Proposition 3.3. There is a smooth family of self adjoint matrices Sp on 
a neighborhood of (p, 0) such that 

i) ReSpP>0, 

(3 ' 15) a) ReSp > K (n+)*n+ - (iip)*np 

This implies the estimates (12.84(1 (12.85(1 which where used in the previous 
section. 
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To analyze H, we use polar coordinates for £ = pQ as in (|2.77|) so that 

(3.16) H(p,C) = P H(p,C,p), H(p,C,p)=H (p,0 + O(p). 

By Lemma l2~T21 for £ G M+ +1 \{0}, # has no eigenvalue on the imaginary 
axis, hence the number N~ of eigenvalues of H in {Hep < 0} is constant. 

We fix a point £ £ £+, that is £ = (r, 77, 7) in the unit sphere with 7 > 0. 
The goal is to construct smooth symmetrizers for H, for (p, p) close to 
(p, f, 0). For convenience we introduce the following terminology. 

Definition 3.4. A smooth symmetrizer for H on a neighborhood uj o/(p, £,0) 
is a smooth self adjoint matrix T> H (p, £,p) such that 

(3.17) Ret H H = J2Vk^kV k , 

where the Vk and are smooth matrices on u of appropriate dimension so 
that the products make sense, satisfying 

i) V k *Vk is definite positive, 

ii) either is definite positive or = 7X^,1 + pX/t^ with Xj; 1 and 
Sfc,2 definite positive. 

Definition 3.5. A family of smooth symmetrizers X K on neigborhoods uj k 
of (p, 0) is called a K-family of symmetrizers for H if there are a decom- 
position 

(3.18) < C N = E H ®E+ 

with dimE~ = iV_ and m{n) — > +00 as k, — > +00 such that for all k 

(3.19) s K (p,c,o) > m( K )n;n+ - it n_. 

where Yi± are the projectors associated to the splitting (|3.18j) . 

Remark 3.6. Recall from MeZu3 that if there is K-family of symmetriz- 
ers, then E~ is the limit of the negative spaces E~(p,(,p) as (p, C,p) tends 
to (p, (, 0) with p > 0. Thus E~ is uniquely determined. On the other hand, 
E + is arbitrary, provided that the the splitting 1)3.18(1 holds: if (|3.19|) holds 
for some choice of E + , then it also holds for another choice for a multiple of 
X K with some other function m(/«). 

We can now state the main result of this paper, which extends |MeZu2j . 
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Theorem 3.7. Suppose that the Assumptions of Section 2.1 are satisfied. 
Assume further that one of the following two condition is satisfied: 

i) all the real characteristic roots (p,r,£) with |£| = 1 satisfy the con- 
dition (BS) of Definition s.^ 

ii) the system is symmetric dissipative in the sense of Definition \2.& 
and the real characteristic roots (p, t, £) with |£| = 1 are either totally 
nonglancing in the sense of Definition \4-3[ or satisfy the condition (BS) of 
Definition \4-V\ 

Then, for all Q G S + , there exists K-families of smooth symmetrizers for 
H(p,C,p) near (p,(,0). 

The condition (BS) ensures that a suitable generalized block structure 
condition, which makes the construction of symmetrizers given in jMeZulj 
work, is satisfied. For hyperbolic problems, it is shown in |MeZu2| that the 
block structure condition is satisfied if the system is smoothly diagonaliz- 
able. In the viscous case, things are more subtle and the generalized block 
structure condition is discussed in details in Section 3. We just point out 
here the following example. 

Theorem 3.8. If (p, f, £) is a semi-simple characteristic root of constant 
multiplicity, then the condition (BS) of Definition \4~Q is satisfied at that 
point. 

Together with Theorem l3.71 this implies Theorem ll.il Finally, we quote 
that the existence of K-families implies the validity of the maximal esti- 
mates when the boundary conditions satisfy the uniform spectral stability 
conditions. 

Theorem 3.9. Suppose that there exists a K-families of symmetrizers for 
H near (p, £, 0) and suppose that the boundary conditions are such that the 
uniform spectral stability condition is satisfied for low frequencies. Then the 
uniform stability estimates ([2.53ft are satisfied. 

Similarly, if the reduced boundary conditions satisfy the reduced uniform 
stability condition then the uniform estimates 1)2 .82j) and 1)2.83)1 hold true. 

3.3 Block reductions 

The advantage of the notion of K-families is that it is independent of the 
boundary conditions. Therefore, their construction depend only on an anal- 
ysis of H. In particular, we can use spectral block decompositions of H. 
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Fix £ G S + . Consider the distinct eigenvalues // of Hq(j>, (). For (p, f , p) 
in a neighborhood of (p, £, 0), there is a smooth block reduction 

(3.20) V~ X HV = diag(^) 

where the have their spectrum in small discs centered at fj, that are 
pairwise disjoints. Equivalently, there is a smooth decomposition 

(3.21) C N = ($E k (p,(,p) 

k 

in invariant spaces for H(pX,p) and .fffc is the restriction of H to E/%. We 
denote by iVj. the dimension of Ej,, that is the size of H^. 

The K-families of symmetrizers are constructed for each block sep- 
arately. If is a K-family for i^fc, it is clear that S K = F*diag(££)V has 
the form l|3.17j) and is a K-family for H. 

When the mode is elliptic, that is when Re fj, ^ 0, the construction of 
symmetrizers is easy (see e.g. |Kre| IUnPi| IMeZulj ). 

Proposition 3.10. Suppose that fj, is an eigenvalue of H (p, Q with Re p k ^ 
0. Then is a smooth family of self adjoint matrices ££ on a neighborhood of 
(p, 0) such that 

i) Re(££# fe )>0, 

(3.22) ii) ReS^>Kld if Re^, fc > 0, 

Re££ > -Id if Rep k < 0. 

Therefore we now restrict our attention to a nonelliptic mode: 

(3.23) Bk = ii d , 

By definition of Hq, this implies that — f + vy is an eigenvalue A of A(p, £) 
with £ = (t/,£ ). In particular, by hyperbolicity, this can only happen when 
7 = 0. By Lemma 12.121 has no eigenvalues on the imaginary axis when 
p > 0, thus the number of eigenvalues in {Re// < 0} is constant. We call it 
Nr. The next definition reformulates Definitions 13 . 41 and 13 . 51 for nonelliptic 
blocks Hk- 

Definition 3.11. A smooth symmetrizer for a nonelliptic block on a 
neighborhood to of (p,(,0) is a smooth self adjoint matrix T,(p,(,p) such 
that, for some C, c > 0, there holds for all (p, C, p) £ oj, 

(3.24) ReE^ fc = 7^1 + ^2, 
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with C, 0) and ^(p, (,0) definite positive. 

A family of smooth symmetrizers on neigborhoods lo k of (p, 0) is 
called a K-family of symmetrizers for if there are a decomposition 

(3.25) E fc (p,C,0) = E fe ©E+ 

wraia dimE^T equal to N7 and m(n) —* +oo as k — > +oo such that for all k 

(3.26) ££(p,C,o) > m( K )(n+)*n+ - (n-)*n-. 

where Ii^ are £/ie projectors associated to the splitting (|3.25|) . 

Given £ = (f , ??, 0) G and a nonelliptic mode = — f, is an 
eigenvalue of A(p,£) with £ = (??,£■)• Therefore, Theorem 13.71 is an imme- 
diate corollary of Proposition I3.1U1 and the following two theorems. 

Theorem 3.12. Suppose that the system is symmetric dispersive in the 
sense of Definition \2. $ suppose in addition that (p,f, £) is a totally incoming 
or outgoing characteristic root in the sense of Definition \4.$ Then there are 
K-families of symmetrizers for the associated block Ht, with E^ = {0} in 
the outgoing case and E^ = C k in the incoming case. 

Theorem 3.13. If(p,f, £) is a characteristic root which satisfies the gener- 
alized block structure condition of Definition \^.2^ Then there are K-families 
of symmetrizers for the associated block . 

4 The generalized block structure condition 
4.1 Hyperbolic multiple roots 

We first recall from |MeZu2j several notations and definitions concerning the 
characteristic roots of the hyperbolic part L. For simplicity, we suppose, as 
we may, that the coefficient of d% is Aq = Id, so that, with notations (|2,3j) . 
L = L. The characteristic determinant is denoted by 

(4.1) A(p,T,Z):=dek(Tld + A(p,t)). 

Definition 4.1. Consider a root (p, r, £) of A(p, r, £)) = 0, of algebraic 
multiplicity m in t. 

i) (p,r, £) is algebraically regular, if on a neighborhood u of (p, £) there 
are m smooth real functions \j(p,£,), analytic in £, such that \j(p,S,) = ~Z 
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and for (p,£) € oj: 

m 

(4.2) ^(P,r,0 = e(p,r,OH{r + X^O) 

where e is a polynomial in r with smooth coefficients such that e(p, t, £) 7^ 0. 

iij (p, r, £) is geometrically regular if in addition there are m smooth 
functions ej(p,£) on lu with values in C N , analytic in such that 

(4.3) A(p,Oej(p,0 = XjMejfat), 

and the e±, . . . ,e m are linearly independent. 

Hi) (p,z,£) is semi-simple with constant multiplicity if all the Xj's are 
equal. 

Case Hi) occurs when X(p, £) is a continuous semi-simple eigenvalue of 
A(p,£) with constant multiplicity near (p, £), such r + X(p, £) = 0. This 
implies that A is smooth and analytic in £ as well as the eigenspace ker(^4 — A). 
In this case, one can choose for {ej} any smooth basis of of this eigenspace. 

If all the roots at (p, £) are geometrically regular, then, locally near (p, £), 
A(p, £) is smoothly diagonalizable, meaning that it has a smooth basis of 
eigenvectors. 

Example 4.2. For the inviscid MHD, the multiple eigenvalues are alge- 
braically regular, but some are not geometrically regular (see jMeZu2j and 
Section |H1 below) . 

The second notion which plays an important role in the analysis of hy- 
perbolic boundary value problems is the notion of glancing modes . Recall 
from McZu2 the following definition. If r is a root of multiplicity m of the 
polynomial A(p, •,£), then by hyperbolicity, the Taylor expansion of A at 
{p-i Hi a ^ the order m — 1 vanishes so that 

(4.4) A(p, t + r, £ + = A m (r, £) + 0(|r, £| m+1 ) 

and A m is homogeneous of degree m. Moreover, A m is hyperbolic in the time 
direction. Indeed, any direction of hyperbolicity for A(p, •) is a direction of 
hyperbolicity for A m . Denote by r + the open convex cone of hyperbolic 
directions fot A m which contains dt. 

Definition 4.3. The root (p,T, £) of A, of multiplicity m, is said nonglanc- 
ing when the boundary is noncharacteristic for A. 

It is totally incoming [resp. outgoing] when the inward [resp. outward] 
conormal to the boundary belongs to T + . It is totally nonglancing if is either 
totally incoming or totally outgoing. 
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Example 4.4. This definition agrees with the usual one for simple roots, 
given by r + X(p, £) = 0. In this case dt + V^A • d x is is the Hamiltonian 
transport field for the propagation of singularities or oscillations and the 
glancing condition d^ d X = precisely means that the field is tangent to the 
boundary. More generally, if the root (p, r, £) of A is algebraically regular, 
then, with notations as in (|4,2|) 

m 

(4.5) A m (r, = e (p,T,|) J] (r + £ Vf A,- <£,£)) 

i=i 

The mode is nonglancing if none of the tangential speed d^ d Xj(p,^) vanish. 
It is totally incoming [resp. outgoing] if they all are positive [resp. negative]. 
In particular, in the constant multiplicity case, all the Xj are equal and they 
are all glancing, incoming or outgoing at the same time. 

In the study of boundary value problems, the dichotomy incoming vs out- 
going plays a crucial role: for instance, for transport equations one boundary 
condition is needed in the first case and none in the second. Using sym- 
metrizers to prove energy estimates, they are constructed in opposite ways. 
The general Kreiss construction also reflects this dichotomy. Introduce the 
following definition: 

Definition 4.5. Suppose that (p,T,£) is an algebraically regular root of A. 
With notations as in Q4.2JI , denote by Vj the order of £ d is a root of order 
of the equation r + Xj(p, £ , . . . , £, , •) = 0, that is the positive integer such 
that 

(4.6) dl d Xj(p,C) = ioia<Vj and := —^XjfaQ / 0. 

We say that Xj is of type I when either Vj is even or Vj is odd and (3j > 0. 
It is of type O when Uj is odd and (3j < 0. 

We denote by Jo [resp. Ji] the set of indices j of the corresponding type. 

Remark 4.6. When (p, f, £) is non glancing, then the all the Vj are equal to 
1, and being of type / [resp. type O] means to be incoming [resp. outgoing]. 
They are all of the same type exactly when the mode is totally nonglancing. 

Remark 4.7. The details of the construction of Kreiss' symmetrizers de- 
pend strongly on being of type I or O, see |Kre| IChFH IMe3j and Section 5. 
There are no reason other than technical why even roots are of type I rather 
than O. 
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4.2 The decoupling condition 

The spectral properties of A(£) are modified by the perturbation B. In 
particular, since the construction of symmetrizers depends deeply on the 
property of being incoming/outgoing, it is very important that the pertur- 
bation respects the decoupling between the different type of modes. 

Definition 4.8. Suppose that (p,r,£) is a geometrically regular root of A 
of order m. Consider a basis {ej} as in (|4.3[) and dual left eigenvectors Ia 
such that 

(4.7) £j (rid + A(p, 0) = 0, l_ d ■ e f (p, £) = 6j, f . 
Consider the and the m x m matrix with entries 

(4.8) 4i'=^ B fe£) e i'fei) 

i) We say that the decoupling condition is satisfied if 

(4.9) B h> = ° when U,f)£(JoXJi)U(JiXJ ) 

where Jq and Ji are introduced in Definition ^. 5\ 

ii) We say that the basis {ej} is adapted to B if 

(4.10) Re5 s >0. 

Definition 4.9. We say that the root (p,r, £) of A satisfies the condition 
(BS) if is is geometrically regular root, satisfies the decoupling condition 
()4.9I) and there is an eigenbasis basis {ej} adapted to B. 

We give several examples and counterexamples. The next result rephrases 
Theorem I3~K1 

Proposition 4.10. (Constant multiplicity) Suppose that (p,T,£) is a semi- 
simple characteristic root with constant multiplicity of A. Then the condi- 
tion (BS) is satisfied. 

Proof. For semi-simple characteristic root A with constant multiplicity either 
Jo or J j is empty so that the decoupling condition (|4.9|) is trivially satisfeid. 
Moreover, it is proved in |MeZulj that (HI) implies that the spectrum of 
B^ is located in {Kez > 0}. Thus there is a basis {ej} in ker(A(p, £) +rld) 
such that HeB^ is definite positive. Next, since any smooth basis {ej} in 
ker(A — A) satisfies (|4.3j) . one can choose it such that ej(p,£) = e_j. □ 
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Proposition 4.11. (Artificial viscosity) Suppose that (p,T,£) is geomet- 
rically regular for %A + B in the sense that there are m smooth functions 
Xj(p,£,p) and m linearly independent smooth vectors ej(p,£,p) on a neigh- 
borhood of (p,£,p), analytic in £, such that \j(p, £, 0) = — r for all j and 

(4.11) {iA(p,£)+pB(p,0)e 3 (p,!i,p) = iXj{p,(, p)e j (p,^ p). 

Then, the decoupling condition is satisfied and the basis {ej|p = o} is adapted 
to B. 

Proof. Alternately, differentiating Q4.3JI with respect to p and multiplying on 
the left by £j/, implies that • = when j ^ j'. Moreover, (HI) implies 
that B\j > 0. □ 

For example, if (p,r, £) is geometrically regular for A in the sense of 
Definition 14,11 and if B = A x Id is an artificial viscosity, then [p, r , £) is 
geometrically regular for iA + B. However, this condition is too restrictive 
for applications, in particular when A and B do not commute. 

Example 4.12. If the root is totally nonglancing, then the decoupling con- 
dition is trivially satisfied since either Jj or Jo is empty. This applies to 
fast shocks in MHD. 

Counter example 4.13. Slow shocks in MHD do not satisfy the decoupling 
condition, see Section |HJ 

The decoupling condition is crucial in the construction of symmetrizers. 
The second condition (14.10(1 is more technical. One could expect that with 
the positivity Assumption (HI), one could always find an adapted basis. 
This is not clear, except for mutliplicity 2 or symmetric systems. 

Proposition 4.14. Suppose that (p, r, £) is geometrically regular of multi- 
plicity m. Assume that either m = 2 or that the symmetry assumption (HI') 
is satisfied. There is a basis {ej} adapted to B. 

If in addition all the eigenvalues Xj are of the same type O or I, then 
the condition (BS) is satisfied. 

The proof is given Section H3 
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4.3 The hyperbolic block structure condition 



We turn back to the construction of symmetrizers for nonelliptic blocks Hj, in 
the splitting (13.20(1 . The construction of K-families is performed in |MeZul| 
provided that can be put in a suitable normal form. This is the so called 
block structure condition. We first review this condition in the hyperbolic 
case, and next extend it to the hyperbolic-parabolic case. 

Consider p and a frequency £ = (r, f), 0) 7^ and a purely imaginary 
eigenvalue (|3.23|) = i£ d of Hq(p,(). Let £ = (77, Then (p,f, £) is a 
root of A. We consider the block H k associated to and denote by E^ 
the corresponding invariant space of H. We use the notations Hk,o(pX) = 
H k (p, C, 0) and Ek.ofo C) = Efcfe C, 0)- 

Definition 4.15. i?fc 5 o the block structure property near (p,Q if there 
exists a smooth invertible matrix V k> o on a neighborhood of that point such 
that V^QHkfiVkfi is block diagonal, 



(4.12) 



V k o HkflVkfi 



Qi 
'•■ 

Q m > 



with diagonal blocks Qj of size Vj x Vj such that : 

QjiPiC) has purely imaginary coefficients when 7 = 0, 



(4.13) 



Qj(P,0 = nJd + i 



1 

'•■ 

'•. '•■ 1 

... 



and the real part of the lower left hand corner of d^Qj(p,Q, denoted by q^, 
does not vanish. 

When Uj = 1, Qj(p,() is a scalar. In this case, (|4.13|) has to be under- 
stood as Qj(p,C) = Mi,; with no Jordan's block. The lower left hand corner 
of the matrix is Qj itself and the condition reads := d^Qj(p, Q 7^ 0. 



Proposition 4.16 (|MeZu2|). If the root (p,f,£) is geometrically regular 
in the sense of Definition ^. 1[ the corresponding block Hk,o satisfies the block 
structure condition. 

Conversely, if Hkfi satisfies the block structure condition with matrices 
V that are real analytic in then the root (p, f, £) is geometrically regular. 
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Remark 4.17. There is a slight discrepancy here between the necessary 
and the sufficient condition, due to analyticity conditions. Definition 14.11 
requires analyticity in £. This is used in the proof of sufficiency. In addition, 
it implies that the block structure condition holds with matrices V that are 
real analytic in £. Thus, there is an "if and only if" theorem. However, 
for the construction of symmetrizers, analyticity of Vk is not needed, this is 
why we do not insist on it in the definition above. In addition, note that 
for fixed p, the existence of C°° eigenvalues and eigenvectors for A, implies 
that these eigenvalues are real analytic in £ and that one can choose analytic 
eigenvectors (see e.g [Sh.il IMal] . The question is to control the domain of 
analyticity as p varies. In applications, for this problem, proving analyticity 
is not harder than proving the C°° smoothness. 

To prepare the hyperbolic-parabolic analysis, we have to review the 
proof of Proposition 14.161 In particular, we reformulate the conditions 
of Definition 14.151 in a more intrinsic way. The choice of a smooth ma- 
trix Vfc 5 o is equivalent to the choice of a smooth basis of E^q, denoted by 
{<Pj,a(p,()}l<j<m',i<a<v r The property (|4.13j) reads 

(4.14) (H o (p,Q- lik )ip jjl (p,() = 0, 

(4.15) (H (p,() - H k )¥j,a{PiC) = i( Pj,a-i(p,,(), 2 < a < Uj. 
With (|3.21|) . there is a unique smooth dual basis ipj^aipX) such that 

^ a ■ K o = 0, 

(4.16) h M 

i>j,a ■ fj',a' = ^j,j'^a,a'- 

Here, K' k denotes the invariant space of Ho(p, () such that = Efc 5 o©lEfc ti- 
lt is the sum of invariant subspaces associated to eigenvalues n k , ^ 

In the basis (fj, a , the entries of the matrix VCq -£ffc,o^4,o are ^j^Hofj^a'- 
The diagonal block structure means that 

(4.17) Vj,a-H'oV 9 i',a' = when j / f. 
The other conditions read: 

(4.18) Re (ip j>a H (p j>a ,) = when 7 = 0, 

(4.19) Re^(^,,^o^,i)(^C) ¥= 0. 

We first show how to compute this quantity in terms of A only. 
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Lemma 4.18. Suppose that Hk,o has a block diagonal decomposition (j4,12l) 
in a smooth basis <f>j jU o/Efc(p, £,0) which satisfies (|4,14|) (|4.15f) . Let 
denote a dual basis satisfying (|4.16|) , The lower left hand corner entry of 
d*fQj(p,Q is equal to the lower left hand corner entry of —idfQj(p,C) and 
equal to 

(4-20) q. = -VWM) ^d\p) m(l>& 

Proof. Let H_ = Hq(j),Q. Then H_ — n k is invertible on K' kQ (p, (). With 
(j4~Ti)) P?T3|) . this implies that 

(4.21) range (H - M fe H) = {VVife C), ■ ■ ■ , V>m> m , (p, C)} -1 , 

(4.22) ker^-^H) = {<p 1}1 (p, () , . . . ,(p m ',i(p,Q}. 

In particular, 

(4.23) (H - ^Id) ^-i = and ^ (ff - ^Id) = 0. 
The entry in consideration is 

Qj(p, = 4>j,v 3 H tp j:1 = ij} jtV . (H - fjLjId)<p jt i + /£ fc <5^,i. 
Therefore, differentiating in 7 and f and using (|2.27f) . implies that 

(4.24) 8^q 3 (p, Q = -id fgj (p, C) = q. 

is given by (|QU|) . □ 

We now discuss how much flexibility there is in the choice of the basis 
ifj^a- Recall that we are considering a purely imaginary eigenvalue fi k = 
of Hq(p, (), so that — r is an eigenvalue A of A(p, £) with £ = (77, £ ). 

Lemma 4.19. Suppose that H^ q has the block structure property near (p, () 
in a smooth basis (fj,a 

and denote by tpj ja the dual basis (|4,16|) . Then, 

i) A is a semi-simple eigenvalue of A(p,£) with multiplicity m equal to 
the number m! of blocks Qj, 

ii) on a neighborhood of (p,(,), there are m smooth eigenvalues Xj(p,C) 
of A(p,£) andm smooth linearly independent eigenvectors ej(p,£), such that 

(4.25) A i b,|) = A, 

(4.26) A(p, i)e 3 (p, i) = Xj(p, |) ei (p, £), 

(4.27) e j (p,i) = ( PjtL ( P X), 
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in) the order of £, as a root of f + \j(p,rj, •) = is equal to Vj, 

iv) denoting by the left eigenvector dual basis of {ej} as in (|4.7|) . 
there holds 

(4-28) l j A d (p)=P J ^(p,Q. 

with Pj := -^-\d v ^ d Xj(p, |) as in (j4.6|) . 

v) the lower left hand corner entry of d^Qj(p,() is 
(4.29) q = -1/Pj G K. 

Proof, a) Define tpj jV j = fj,vj an d for a < Vj 

(4-30) £/,a(Pi C) = -i{H (p, C) - (jQv>j,i/j- 

By dH^mUD, there holds 

(4.31) £j>(p>C) = <Pj,a(p,£)- 

Moreover, in the new basis <^j, a , the matrix of Qj has the form 



(4.32) Qj = i£jd + i 



Thanks to (|4,31|) . the dual basis {tpj, a } associated to {</5j, a } also satisfies 
^jaiPiC) = ^jaiPiC)- This implies that the lower left hand corner of 
dx/Qj(p, Q is unchanged in the new basis. 

b) Consider the determinant 

Aj(p, C, id) = det (&Id + iQj-(p, C)) • 

It is independent of the basis {ipj, a } or {ipj a}. Thus, it is real when 7 = 
and vanishes at (p,C,€ d )- Moreover, Q4.13JI implies that 

As a byproduct, using also (|4.24|) this shows that 

(4.33) q.£R thus q . = Re q . = q\ \± 0. 
—j —j —j J 

In particular, the implicit function theorem implies that there is a smooth 

function Xj(p, £), in a real neighborhood of (p, £), such that Aj(p, £) = —f 

and for £ = (f, 77, 0): 

(4.34) Aj (p, C, | d ) = ay (p, C, 6) (f + A i (p, I)) 
with aj{pX,i d ) / 0. 



/* 


1 . 


. 0\ 
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c) Consider next the eigenvector equation 

(4.35) (i d ld + iQ j (p,C)h=0. 

By 1)4.32)) . in the basis {ipj, a }, the Vj — 1 first equation determine the last 
Vj — 1 components of ej 

(4.36) (ej-) a = (Id " I d ) a_1 (ei)i, a > 2. 

Substituting these values, the last equation is a scalar equation equivalent 
to Aj = 0. Introduce 

Cj (Pi »7> I) = ( - A j (p> £) > *7> °) , 

and 

(4.37) ej (p, i) = ^,i(p, C) + - ^'"Vi.aCp, 0- 

a=2 

This vector is smooth and satisfies 1)4.35)1 . thus 

(A(p,£) - Xjfaftltfejfai) = A d (p){iH G {p,Q) + &Id) ei (p,£) = 0. 

Moreover, the ej(p, £) = y_j i(p, C) are linearly independent. 

d) Bv fOP) , for ( = (f,7/,0), there holds 

det (fid + A(p,i)) = det(A d ) det (i# (p, C) + &H) 

m' 

= «(p,r, |) JJ(f + \j(p,i)) 
i=i 

where a(p,r, £) 7^ and m' is the number of blocks Qj. This shows that — f 
is an eigenvalue of algebraic order m' of A(p,£f\. By step c), the geometric 
multiplicity is at least m' , implying that —f is semi-simple of order m! . 
Moreover, by 1)4.14)1 . there holds 

showing that £ is a root of multiplicity Uj of Aj, thus of r + Xj(p, fj, £) = 0. 
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e) Let ij satisfy (|£7|) . Thus 

Range(H (p, Q - /f fc Id) = A^ip) Range (fid + A(p,[)) 

= Aa 1 (p){£ 1 ,...,lm} ± - 
Comparing with (|4.21f) . this implies that 

(4.38) span{tpj jUj (p, Q, 1 < j < m} = spanj^, 1 < j < m}. 

For a £ {1, ... , Uj}, introduce 

(4-39) %|B = _i_fl« r ie,( g ,£). 

Because £ d is a root of order Uj of f + Xj(p,f/,^) = 0, the definition (j4,37|) 
implies that 

%,a = £?,a(P,C) = ¥j,a(P,0 for 1 < a < Uj. 
In particular, (|4.16|) implies that 

(4.40) tpj,^, (p, C) • e jtV . = il)jt tVjl (p, C) • (p, C) = $j,j>- 
Differentiating the equation 

(4.41) (A(p, i) - Xj (p, 0) ej (p, = 
with respect to ^ and at order Uj yields 

(fid + A(p,£))d£ej(p,i) = -VjAaipjd^ejfai) + dg'A^l^f). 

Multiplying on the left by iji annihilates the left hand side, implying 

£j,A d (p)ej iVj (p,() = 0j£j, ■ ej(p,i) = (3j5 fij . 

By (|4.38|) . the i_jA d and tpj,^ span the same space. , Therefore, comparing 
with (|4~471 implies that ijiA d (p) = (3jipj, tV .,(p,(). 

f) By (ET2U]) and (I4~2%1) . we have 

-PjQj = lj<Pj,l(P'i) = lj e j(P>t) = L 
The proof of the lemma is complete. □ 
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Remark 4.20. This lemma is a variation on the necessary part in Propo- 
sition 14,161 (see McZu2 ), with useful additional remarks. It shows that 
the block structure condition is closely related to a smooth diagonalisa- 
tion of A. Conversely, if one starts from a smooth basis e» and a root of 
f+Xj(p, £) with (|4.6|) . one constructs a basis ifj ja such that tpj. a (p, C) is given 
by ()4.39|) , using an holomorphic extension of ej to complex values of ^ (see 
MeZu2 ) . Lemma 14.191 implies that the change of bases which preserve the 
block structure form are linked to change of bases which preserve the smooth 
diagonalization of A. 

The construction of K-families of symmetrizers for the blocks Qj is per- 
formed in [Krel |MajH lMe3j . The sign of Pj and the parity of Vj play an 
important role. Hyperbolicity implies that Hq and thus the H k and Qj have 
no purely imaginary eigenvalues when 7 > 0. Denote by Eg . the invariant 
space of Qj associated to the spectrum in {Re u < 0} since the definition of 
the limiting space Eg . Recall that the limit space at (p, Q is 

(4.42) Eg\ = &i x {Op -1 '.* 

with 

{Vj/2 when Vj is even, 

(za, + 1)/2 when Uj is odd and > 0, 
(yj — l)/2 when Uj is odd and (3j < 0. 

Remark 4.21. As a corollary, we have the following characterization of the 
sets Jo and Jj: 

( j E J 1 if fj is even or Vj is odd and q^j < 0, 

[ j £ Jq if Vj is odd and > 0. 



(4.44) 



4.4 The hyperbolic-parabolic case 

We still consider a block associated to a purely imaginary eigenvalue 
()3.23j) . In the next section, we show that the following technical conditions 
are the natural one for the construction of Kreiss symmetrizers. 

Definition 4.22. has the generalized block structure property near (p, 0) 
if there exists a smooth invertible matrix on a neighborhood of that point 
such that 

( Bn ■ ■ ■ Bim 



(4.45) Vr l H k V k 



( Qi ■ ■ ■ 
\ ••• Q„ 



+ P 



\ B m \ • • • B 7l 
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where the Qj(p,C) satisfy the properties of Definition \4- i3| Moreover, the 
m x m matrix B 9 with entries B\ ■, equal to the lower left hand corner of 

Bj,j'(Pi()fy satisfies 

(4.46) B ii> = ® when (j, f) e (J x J/) U ( J 7 x J Q ) 

where Jq and Jj are defined by (|4.44() and there is a real diagonal matrix 
, with entries cOj such that 

(4.47) d)q) >0, Re£> b 5 b >0. 

We show that these conditions are related to the condition (BS) of Def- 
inition formulated on the original system. We need first a more detailed 
form of the block reduction H in ()2.28|) . Introduce the following notations: 

d-i 

(4.48) B^(p,() ■= VjVkB jtk (p), 

j=i,k 

d-l 

(4.49) B^(p,C) :=J2vo{ B j,d(p)+ B d,j(p)) 

Lemma 4.23. One can choose the matrix V in (|2,28|) such that there holds 

(4.50) H(p, C) = H (p, C) - H x {p, + 0(|C| 3 ) 
where 

(4.51) H x = Af (X* - iB* >d H - B d4 Hl) , 

Proof. Direct computations show that the kernel of G(p, 0) is x {0} and, 
using that A d is invertible, that kevG(p, 0) nrangeG(p, 0) = {0} This shows 
that is a semi-simple eigenvalue of G(p, 0). 

If fj, is a purely imaginary eigenvalue of G(p, 0), then is an eigenvalue 
of iA(p, £) + B(p,(;) with £ = (0,—ifi). By Assumption (HI) this requires 
that £ = 0, thus [i = 0. This shows that the nonvanishing eigenvalues of 
G(p, 0) are not on the imaginary axis. 

This implies that there is a smooth matrix V(p, () on a neighborhood of 
(p, 0) such that ([2.28ft holds with H(p, 0) = and P(p, 0) invertible with no 
eigenvalue on the imaginary axis. 
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The image of the first N columns of V is the invariant space of G, and 
H is the restriction of G to that space. At £ = this space is kerG, and 
performing a smooth change of basis in C^, we can always assume that the 
first N columns of V are of the form 



(4.52) VjipX) 



( Idjvxiv 

\w(p,o 



with W of size N' x N vanishing at £ = 0. This implies (|2.3U|) . 
By (H2j GVi = VjH, hence MVj = G d ViH and 

M = -AH + B d WH, W = JR. 

Therefore, 

(4.53) M = -AH + B d JH 2 = -AH + B d4 H 2 . 

Taking the first order term at £ = shows that the first order term in Hq 
in H satisfies 

d-l 

(it + 7 )id + l m A 3 = -Mp)h 

and hence is given by (|2.27|) . The second order term Hi in H satisfies 

i?*,* = —A d H\ + iB* d Ho + B d ^ d JH§ 
implying (|43Uj) and (|43Tj) . □ 
Parallel to Lemma 14.181 we can now state: 



Lemma 4.24. Suppose that the matrix of H^ is given by the right hand 
side of (|4.45|) in a smooth basis ipj ta o/E^(p, p) which satisfies (|4,14|) and 
(|4.15() /or p = 0. Zei denote the dual basis of {ej = <fj t i} satisfying 
(|4.7jl . TTie entries of B 9 are 

(4-54) B^, = -1^.^,1)^(^,0). 

Proof. In the block reduction (|4.45l) , the lower left hand corner entry of the 
(j,f )-block is 
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Differentiating in p and using the relations (|4,23|) yields 



-b)# = d p hjAp,c,o) = -± jn B(pX)<p j}1 , 

where ipj u , and (p. 1 stand for the evaluation at (p, C, 0) of the corresponding 
function. Using the explicit form of B and the relations 

we obtain 

With (|Qg|) . this implies l|131|) . □ 

Theorem 4.25. If(p,f,£) is a geometrically regular characteristic root of 
A which satisfies the condition (BS) of Definition \4 ■ !\ Then the associated 
block Hk satisfies the generalized block structure condition. 

Proof. Since (p, f, £) is geometrically regular, the hyperbolic part H^fl sat- 
isfies the block structure condition. Moreover, if e,- is a basis analytic in 
£, there is a basis <pj ja such that ipj !a (p,() = ej(p,S t ) (see Remark |4,20l or 
|MeZu2p . By Lemma H~2ll (EP1) is equivalent to EST)!! . 

If once can choose the base {ej} such that (J4.1U|) holds, then choose 
= -Pj and by gjjH) and (jHSl|) there holds djgj = 1 so that L>5 b = £?# 
satisfies (|4~IT|) . □ 

Remark 4.26. Conversely, if the generalized block structure condition 
holds with matrices which are real analytic in then, by Prop osition 14 . 1 6l 
(P)Z)0 is geometrically regular. By (j4.54|) . (|4.4(i|) is equivalent to the de- 
coupling condition (|4.9|) . Moreover, ()4.47j) implies that there is a diagonal 
matrix with positive entries d!- = d^/q^- such that Re D^B$ > 0. Consider the 
diagonal matrix C = = diag(cj) and the new basis e~j = CjCj. The 

new dual basis is lj = cj 1 ^ and the new matrix & is C~ l B^C = CD^B^C 
and therefore Ke& = CRe (D^B^)C is definite positive. 
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5 Symmetrizers 



In this section, we prove Theorems 13 . 1 21 and 13 . 1 31 We are given a frequencies 
C = (r^jO) an d a purely imaginary eigenvalue p^ = — i£ d of Ho(p,Q, so 
that (p,f,£, with £ = (jjt^j) is a ro °t of the characteristic determinant A, 
of multiplicity m. Our goal is to construct K-families of symmetrizers for 
the block Hk{p,C,-,p) associated to pk- 



5.1 Proof of Theorem I5T51 

We assume here that (p, f, £ is geometrically regular and satisfies the condi- 
tion (BS). 

We follow closely |MeZulj (Lemma 4.11 and Appendix A therein. See 
also |Me3j ) where the constant multiplicity case is studied. In this case, 
all the blocks Qj are equal and thus have the same dimensions but more 
importantly, all the eigenvalues are of the same type O or /. So we review the 
main steps of the construction and indicate where the proof of [MeZull IMe3| 
has to be modified. 

In the block reduction ()4.45|) of Hk, we choose the symmetrizers to 
be block diagonal: 



(5.1) 



h 







••• S£ 

E? + E?(p,{)- 



/- /•]' - ipFf 



where Ej and E* are real symmetric matrices, and F* and F* are real and 
skew symmetric. Moreover, E? , F? and F? are constant, E* depends only 
on (p,f,fj) and the Ef have the special form 



Ef 











and Ef(p,() = 0. 

The block structure condition implies that 



(5.2) 



H k = diag(Q i | 7=0 ) + 7diag(5 7 Q j | 7=0 ) + pB\ p= qO{^ 2 + p 2 
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is a symmetrizer for H^, on a neighborhood (depending on k) of (p, (, 0), 

if 

(5.3) Re + E«)Q jh=0 ) = 0, 

(5.4) Re (Ep^QjipX, 0) - iF?Qj(p,Q) > 0, 

(5.5) Re (diag(E;)£(p, C, 0) - idiag(F/Q,(p, £))) > 0. 



Moreover, the condition (|3,26|) reads 

|2 Itt-„..|2 



(5.6) (£^w, w) > Ci ( As|nj«?|' a - |njt« 

where 11^ is the projection onto E^ in the decomposition C Uj = E~ © E^", 
where 

(5.7) Ej = C'i x , E+ = {0}^ x , 

with v'- given by ()4.43|) . 

Before starting the construction, we note that Ralston's lemma |Ra| (see 
also |MeZuH iMeS] ) implies that one can perform an additional change of 
basis Id + pV such that the matrices Bjj> in (|4.45j) are of the form 



(5.8) B jtf (p,() 



( * 
: 



This does not affect the previous choices, made at p = 0. Next, we introduce 
nome notations. A vector w S C Nk = ©C^', is broken into m blocks Wj G 
C-ij with components denoted by Wj A . We now proceed to the construction 
of the symmetrizers. 

a) One first choose the Ej such that ()5.6|) holds and 

(5.9) Re {l-: , pLQ J [p.Q l r J .u-) > 2\w j:1 \ 2 - C^w'jf 

(5.10) Re (dmg(E*)B(p,(,0) Wj , Wj ) > 2\w* A \ 2 - C' K \w',\ 2 . 

with Wj t i denoting the first component of Wj G C Uj and w'j £ C^' -1 denotes 
the other components and w* t i £ C m is the collection of the first components 
Wj i while denotes the remaining components. 
Note that 

Re {Efd^Qjip^w,™) = e^Kil 2 + 0(|^mi), 
Re(diag{E^)B(p,C,0)w j ,w j ) = Re (EpB^w*,!, w M ) + O(KIH), 
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where E is the m x m diagonal matrix with entries e^. Moreover, the 
decoupling condition (|4.9|) implies that B b has a block diagonal structure: 
ordering the base {ej} according to the type I or O, with obvious notations 
there holds : 

< 5 -"> M B o ! 4) 

Similarly we note E h = diag(£^,.E^) and and are satisfied if 

(5.12) el iq ) > 3, Re£^>3Id, ReE^ B b > 3Id. 

On the other hand, to satisfy (|5.6|) . one chooses the e^ a inductively, 
starting from a = 1, but this choice depends on the type of the eigenvalue. 
Remember also from (|4.29|) that f3j = —1/q^. According to |Kre| IChPij or 
Lemma 8.4.2 in Me3 , the e^ a are chosen as follows. 

1) If Xj is of type / , then = e^i is taken 0(1), independent of k, 
and the e K - a for a > 2 are chosen successively and depend on k. In particular, 
when Vj is even, e* 2 > c/c. When Vj is odd and /3j > 0, then < and 
e^i < 0; when Uj > 3, then 3 > ck. 

2) If Xj is of type O, that is Vj odd and > 0, one chooses > ck 
and the other e^ a are chosen inductively. 

By assumption, there is a diagonal real matrix D 9 = diag(Z)^, Dq) such 
that 

d)q)>Q, ReD b jB b j>0, ReD b o B b o >0. 

Therefore, there is a positive constant c such that if we choose = cdj 
when Aj is of type / and e*i = cntfj when Xj is of type /, the condition 
(t5~T2l is satisfied. Next, according to jKrel IChPil IMe3] . we can choose the 
e K j a for a > 2 such that the inequality (j5.fi j) is also satisfied. 

Remark 5.1. The construction above shows that the conditions of Defi- 
nition 14.221 are more or less necessary for the construction of K- families of 
symmetrizers. First, the different magnitude in k of for different types 
forces the decoupling (j5,ll|h that it condition (|4.9|h Second, a spectral con- 
dition on B^ is not sufficient in general to insure the existence of a diagonal 
matrix E b such that (j5.12j) holds. This indicates that condition (|4.47|) is 
also necessary for the construction above. 

b) Once the matrices Ej are chosen, the construction goes on as in 
|MeZii1llMe3] . We omit the details. By ijOSIl . Re (EjQj(p, ()) = 0. Next, 
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using the implicit function theorem and the property that -?Qk is re& l when 
7 = 0, the real symmetric matrix Ek(p,f,fi) is chosen so that such that 
Re (E j + E J )(Q jh=0 ) = 0. 

Since Fj is real and skew symmetric, there holds Re — iFjQj(p,<^) = 
HeFjJj where Jj is the Jordan matrix in (|4.13f) . One can choose Fj such 
that 

Re(FjJjWj,Wj) > -\w jt i\ 2 + (C+ l)|u^| 2 . 

where C is the constant in (|5.9[) , Adding to (|5.9|) implies (|5.4j) . 

Similarly, Re — iFjQj(p, = HeFjJj and one can choose Fj such that 

Re(FjJ jWj ,Wj) > -\w jtl \ 2 + (C + l)\w'j\ 2 . 
where C is the constant in (|5,1U|) . implying ()5.5|) , 

5.2 Proof of Theorem HTT21 

We now assume that the symmetry hypothesis (HI') is satisfied and that the 
root (p,f,£, is totally non glancing. In MeZu2 , symmetrizers forijfc(p, 0) 
are constructed. We show that they also symmetrize Hk{p, p) when p > 0. 

In |MeZu2j . it is proved that the nonglancing condition implies that the 
multiplicity of pk as an eigenvalue of Hq(pX) = H(p,(,0) is equal to m. 
Denote by the N xm sub-matrix of V which corresponds to the block H^. 
Therefore, for (p, p) close to (p, 0), the corresponding invariant space of 
H h is E fc (p, C, p) = V k (p, C, p)C m "and 

(5.13) V k H k = HV k 

Recall that (p, C,p) is the negative space of Hk for ( 6 S + , p > with 
7 > + p > 0. 

Lemma 5.2. iij If(p,T,£) is totally incoming, then, for (p,C) * n a neigh- 
borhood of (p,C), Efc (p, C)~= C m . 

m) If (p,Z,^) is totally outgoing, then, for (p,() *n a neighborhood of 

(p,q, ^{p,o = {o}. 

Proof. The dimension is constant for 7 > + p > 0, and the result is proved 
in |MeZu2| when p = 0. □ 

By assumption, there is a definite positive matrix S(p) such that the 
SAj are symmetric. 
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Lemma 5.3. The symmetric matrix 

(5.14) £fc,o(P>0 = -V k *(p,(,0)S(p)A d (p)V k (p,(,0) . 

is a symmetrizer for H k on a neighborhood of (p, 0). More precisely, there 
holds 

(5.15) Re^ k H k = 7 i?i + /9i? 2 

toit/i Si(p, C)0) and £2(^5 C, 0) definite positive. 

In addition, S&(p, (, 0) is definite positive [resp. negative] when the mode 
is totally incoming [resp. outgoing]. 

Proof. According to (|.3.16[) . there holds 

H(p,C,p) =H (p,()+pH'(p,C,rho). 

Using (|5. 13|) and the definition (|2.27f) of Hq, one obtains the identity (|5.15|) 
with 

(5.16) Ri = V k *SV k , 

(5.17) R 2 = V k *(ReSA d H')V k . 

Because S is definite positive, i?i also has this property. Next, Lemma I4.2HI 
implies that H'(p, (, 0) = —Hi(p, Q with H\ given by (|4.51|) . Since Hq(j>, () = 
^fcld = —i^ k ld on Efc(p, C, 0), there holds 

H'(p,(,0)V k (p,(,0) = -A d l (p)B{p,l). 

Therefore, at the base point (p,cz, 0), there holds 

R 2 (p,cz,0) = V k *(ReSB)V k . 

The assumption (HI') implies that SB is definite positive on the space 
Efe(p, CjO) = ker(^4(p,£) + rid), implying that i?2 is definite positive at 
(p,cz,0), hence on a neighborhood of that point. 

That C,0) is definite positive [resp. negative] when the mode is 

totally incoming [resp. outgoing] is proved in |MeZu2| . □ 



With Lemma 15.21 this implies that 
(5.18) = 
are K-familes of symmetrizers for H k . 



Sfc in the incoming case, 
nT, k in the outgoing case. 
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6 Further remarks and examples 



6.1 Adapted basis. Proof of Proposition [4.141 

In this section, we always assume that Assumptions (HI) is satisfied. Con- 
sider a geometrically regular root (p,f, £) of A. We show that there are 
eigenbasis {e^} satisfying Q4.3|) which are adapted to B, in the sense of Defi- 
nition |4~%1 either when the multiplicity is 2 or when the system is symmetric. 

1. The case of multiplicity two 

Projecting on the 2-dimensional invariant space of iA(p,£) + pB(p,£) 
associated to the eigenvalues close the i\± and 1X2, we are reduced to consider 
2x2 matrices 

(6.1) iA(p,i) + pB(p,i,p) with A = 

The Assumption (HI) implies that the spectrum iA + pB is contained in 
Re A > cp, for (p, £) close to (p, £) and p S [0, po], for some po > 0. Changing 
£ to — £ and using (HI) near — £, we see that the spectrum ±iA + pB is 
contained in Re A > cp. 

We show that, changing the base {ei,e2} if necessary, one always meet 
condition (|4.1U|) . 

Lemma 6.1. With assumptions as above, there is a smooth change of bases 
preserving (|6.1|) . such that KeB(p,^) is definite positive. 

Proof. The constant multiplicity case Ai = A2 being already treated, we 
assume that Ai 7^ A2 on any neighborhood of (p, £). In this case we are 
limited to consider diagonal change of basis and we prove that there exists 
a diagonal matrix D, such that 

(6.2) Re (DB(p, |, 0)D- 1 ) > 0. 

a) Recall that there is c > such that the spectrum ±iA + pB is 
contained in Re A > cp. We first show that for all t G M the spectrum of 

(6.3) (° ») + B( B |,0)-£ld 

is contained in Re A > 0. If not, there are t, p\ > and a neigborhood to of 
(p, £ ) such that 

(6-4) (° 



U A 2 i' 
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has an eigenvalue in Re A < c/2 when (p, £) G u; and p G [0, pi]. There 
is G u; such that A20',!') - Ai(p',|') = ii / 0. Choose p G [0,pi[ 

such that p\t\ < \ti\. By continuity, since A2 — Ai vanishes at (p, £), there 
is (p, £) G cj such that A2O5, ~~ ^i(p,Q = Therefore the matrix 

±L4(p, £) + pB(p,£,p) has an eigenvalue in {Re A < pc/2}. 

b) Consider a matrix 

f a 6 N 

c dj 

Its spectrum is contained in {Re A > 0}, if and only if 

Re(a + d)>0 and |Re (a//)| 2 < Re (a + df 



where / = (a - d) 2 + 46c. Since |Re( v / 7)| 2 = |(|/| + Re/), the second 
condition reads 

|Im/| 2 < 4(Rea + Red) 2 ((Rea + Red)) 2 - Re/), 



or 



(Re (a - d)lm (a - d) + 21m (6c)) 2 < 

(Re a + Re d) 2 ((Im a - Im df + 4Re aRe d - 4Re (6c)) 



or 



ReaRed((Ima — Imd) 2 — (Re a — Red)Im (6c)(Ima — Imd) 
( °" )) +(Rea + Red) 2 (ReaRed-Re6c) - |im(6c)| 2 ) > 0. 

We apply this criterion to the matrices ()6.4|) . In this case, when t varies in K 
the coefficient Im (a — d) varies from —00 to +00 while the other coefficients 
are fixed. Therefore, if the corresponding inequality (|6.5|) is satisfied for all 
t, then Re a + Re 6 > 0, Re aRe d > and 

(Rea-Red) 2 |im(6c)| 2 < 

4ReaRed((Rea + Red) 2 (Re aRe d - Re 6c) - |Im(6c)| 2 ). 

Thus 

|Im (6c) 1 2 < 4Re aRe d(Re aRe d - Re 6c) 

and 

\bc\ +Re(6c) < 2ReaRed. 
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Denoting by bj^ the entries of B(p,£, 0), we see that the spectral condition 
of step a) implies the following conditions: 

(6.6) Re6n >0, Re6 22 >0, |& 12 &2i| + Re (&12&21) < 2Re&nRe6 22 . 
c) Similarly, we note that the condition HeB(p, £,0) > is equivalent 



to 

(6.7) 



Re 6n > 0, Re 6 22 > 0, 

\h2 1 2 + I&21 1 2 + 2Re (612621) < 4Re6nRe6 22 . 



With D =(l y), the 

conjugation DBD 1 changes B to B' with 

b' n = b 11 , 6' 22 = 6 22 , 6' 21 = (56 2 i, b' 12 = -b 12 . 
For all e > 0, one can choose (5 such that 

I 6 'i2l 2 + l ft 2il 2 < I&12&21I +e 

and therefore, (|6.6j) implies that there is 5 such that (|6.7j) holds for the 
b' jk - □ 

2. Symmetric systems. 

Lemma 6.2. Suppose that (p,f,£.) is geometrically regular and Assumption 
(HI') holds. Then one can choose the eigen-basis {ej} such that ()4.9|) holds. 

Proof. Denote by S the symmetrizer. We show that one can choose the 
eigen-basis ej such that 

(6.8) *ej(p,|) S(p) eji(p,£) = 5 jd >. 

In this case, ij = t ej(p,^)S(p) and 

(6-9) B)., = * ei (p,|)5(p)B(£,|) e/ (p,|) 

showing that Re B* is the restriction of Re (S(p)B(p,£)^ to the space spanned 
by the ej(p, £), and hence positive. 

To prove (|6.8|) . consider the partition of {1, ... , m} into susbets J a such 
that j and f belong to the same class J a if and only if Xj = \'- on a 
neighborhood of (p, £). Denote by ¥ a (p, £) the space spanned by the ej(p,^) 
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for j G J a . Then, near (p,£), A(p,£) = A a Id on this space, where A a is the 
common value of the \j for j £ J a . Thus, locally, one can find a smooth basis 
of ¥ a , analytic in £ and orthonormal for the scalar product S(p). Collecting 
these bases of F a , (|6.8f) holds when j and j' belong to the same class J a . 

When j and j' do not belong to the same class J a , there is a sequence 
(p n ,£, n ) converging to (p, £) such that \j(p n ,C n ) <V'(P n The symme- 
try of S(p n )A(p n , £ n ) implies that 

t e j (p n ,?)S(p n )e j ,(p n ,Z n ) = Q. 

Therefore, passing to the limit, we see that (|6.8j) is also satisfied when j and 
j' do not belong to the same class J a . □ 

6.2 Discontinuity of the negative spaces E 

We show the decoupling condition (|4.9|) is necessary for the continuity of 
E _ (p, C, p) at p = 0. Before stating the result, we make the following remark. 

Lemma 6.3. Suppose that (j),f, £) is geometrically regular and nonglancing. 
With notations as in flOJ, (|4~%)) Zei = d^ d Xj(p,[) ^ 0. T/ien, tftere 
is c > suc/i i/iaf /or all p > and t £ 1, i/ie spectrum of itdi&g((3j) + pB^ 
is contained in {Re /i > cp} . 

Proof. We fix p = p and forget it in the notations. For £ close to £, consider 
the invariant space of iA(^) + pB(£) associated to eigenvalues close to — r. 
In the basis {e^}, its matrix is 

(6.10) idiag(A i (0)+A>B(^/>), 

with .B(£, 0) = The Assumption (HI) implies that the spectrum of this 
matrix lies in {Hep > cp}. 

Adding rid, we can assume, without loss of generality, that Aj(£) = 0. 
Taking t > 0, £ = £ + ±(0, t) and p = icr > the matrix in in (jfi.K)J) is 

tM(i, o") = t( ± idiag(/3,) + <rB tt + 0(t)) 

and the spectrum of M(t,a) lies in {Re/i > ca}. Letting t tend to zero, 
implies that the spectrum of M(0, cr) is also contained in {Re p > ccr} and 
the lemma follows by homogeneity. □ 

Corollary 6.4. If (p, f, q) is geometrically regular and nonglancing, then 
for all 7 > and p > 0, with 7 + p > 0, the matrix diag(/3~ 1 ) (7H + pB^) 
has no eigenvalues on the purely imaginary axis. 
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Consider £ = (f , 77, 0) 7^ and a purely imaginary eigenvalue \i = i^ d of 
Ho(p, ()■ Let £ = (77, £ ). Then r, £) is a root of A. We denote by the 
block associated to /j and, for /? > 0, we denote by E^~ (p, C, p) the negative 
invariant space of 

Proposition 6.5. Suppose that (p,r,£) is geometrically regular and nonglanc- 
ing and suppose that there exist j £ Ji and j' € Jo such that 

(6.11) flj^. / 0. 

T/ten i/ie negative space E7(p, £,p) has no limit as (C,p) — » (C,0). 

In particular, there are no smooth K-families of symmetrizers for 
near (p, Q . 

Proof. By Lemma 14.181 and Lemma 14.241 the block decomposition (|4.45|) 
implies that in a suitable basis 

(6.12) H k (p, f , 2, 7, p) = -diagOS" 1 ) ( 7 Id + + 0( 7 2 + p 2 ) . 
Denote by IE - (7, p) the negative space of fffc(p,T, 77, 7, p). We show that 

(6.13) limE _ (7,0) / limE~(0,p), 

7— >o p— >o 

which implies that IE - (7, p) has no limit as (7, p) — > (0,0). 

Consider first the case where p = 0. Then, ()6.12j) implies that the first 
limit in 1)6. 13j) is the space Ej spanned by the vectors ej of the basis such 
that (3j > 0, that is such that j £ Jj. 

On the other hand, Corollary 16.41 implies that B° = — diag(P~ 1 )B^ has 
no eigenvalue on the imaginary axis. Therefore, the second limit in (|6.13l) is 
the negative space E~ b of B 9 . If it were equal to Ej, this would mean that 
E/ is invariant by B , thus by = — diag(Pj)B° , which contradicts (|6.11j) . 

By |MeZu3j , the existence of smooth iT-families of symmetrizers implies 
that the limit of E^~ at (£, p) exists, and is equal to the space EJT of Defini- 
tion Therefore, (|6.13l) implies that there are no smooth if -families of 
symmetrizers. □ 

6.3 Viscous instabilities 

Consider boundary conditions as in Assumption 12.91 When the negative 
space E~ is not continuous in (C, p), then the Evans function is likely not 



60 



continuous and one can expect that the low-frequency uniform stability con- 
dition for the viscous problem is strictly stronger than the similar condition 
for the inviscid problem. In particular, the inviscid problem can be strongly 
stable while the viscous one is strongly unstable. We illustrate here this 
phenomenon on an explicit example. 

1. An example. Consider the system 



(6.14) 



J (d t + d y )ui + d x u 2 = efj,Am, 
\ (d t + d y )u 2 + d x ui = evAu 2 . 

Taking linear combinations and changing e, the system is equivalent to 

(6.15) (d t + d y )ld + Ad x - eBA, A=(l , B=(l a Y 

with a = \u — y\l(v + /u) G [0, 1[. This system is symmetric and satisfy the 
assumptions (HI) and (H2). 

The hyperbolic part is diagonal: the eigenvalues are 

(6.16) Ai = + A 2 = 77-£. 

They cross on the line ^ = and are trivially geometrically regular since the 
system is already in diagonal form. One of the eigenvalue is incoming, one 
is outgoing. The decoupling condition (J4.9|) is satisfied if and only if a = 0. 
In the sequel, we assume that a > 0. 

2. Boundary conditions. Next, consider boundary conditions for (|6,15|) : 

(6.17) uy x=0 + er^ui^o = 0. 

We first compute the limiting inviscid boundary conditions, using boundary 
layers. The bounded solutions u = w(x/e) of (|6.15j) are 

(6.18) w(z) = u + e zB ~ lA h, h G E B _ 1A , u £ C 2 . 

where K^_ 1A is the negative space of B~ 1 A. Therefore, u is the endpoint of 
a profile which satisfies the boundary condition (|6.17jl . if and only if 

(6.19) u £ (ld + TB~ 1 A)E B _ 1A . 

Note that given any complex number c, one can choose F such that this 
boundary condition reads 

(6.20) ui = cu 2 
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3. Low frequency stability. The first order system (j2.21|) reads 
(6.21) d z U-G(C)U, G(0=( aB J +lfljd B \), 

with C = (r, rj, 7) and a = 7+i(r+r/). Perform the small frequency reduction 
(|2.28f) . using the change of unknows 




Then, by Lemma 14.231 there holds 

^ = (0 p) 

with P(0) = and 

(6.22) H(C) = -a A + (a 2 - r, 2 )AB + 0(|C| 3 ), 

Since V(0) has the triangular form (j2.30j) . we see that the boundary 
condition reads 

(6.23) u H + f(()u P = 0, f(0) =r + A~ 1 B. 

The Evans condition is violated at £ if there is up 6 Ep(£) and rip G Ep(C) 
satisfying this boundary condition. The negative space of P(C), Ep(£) is 
smooth in £ and equal to E^_j^ when £ = 0. Thus, the Evans condition is 
violated at £ if and only if if there is uh £ Ep(£) such that 

E-(C)nf(c)E;(C)^{0} 

Since A~ X B = (B^A)' 1 , there holds 

f(0)Ep(0) = (ld + TB- 1 A)E B _ 1A . 

Comparing with (|6.19l) and (j6.2Uf) . we see that for £ small, the space r(£)Ep(£) 
is generated by t (c((), 1) where c(() is a smooth function such that c(0) = c. 
Therefore, the Evans condition is violated at £ if and only if 

(6.24) G Ep(C). 
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Remark 6.6. Using the terminology of MeZu2j, the analysis above shows 
that the reduced boundary condition for the hyperbolic part H(Q reads 



(6.25) ui = c(C)«2- 

Taking £ = in this equation, we recover that (|6.2()|) is the natural limiting 
boundary condition for the hyperbolic operator Hq . 

Proposition 6.7. There are choices of a and V, such that 

i) the inviscid problem ()6.15j) for e = with the boundary condition 
(|6.2()|) is maximal striclty dissipative thus uniformly stable, 

ii) the viscous problem with boundary conditions (|6.17|) is strongly un- 
stable for small frequencies, in the sense that the Evans functions vanishes 
for arbitrarily small frequencies £ with 7 > 0. 

Proof. The matrix 

(6.26) 5 = ( I ° ) , s > 



is a symmetrizer for the inviscid problem. If 

(6.27) |c| 2 < s, 

the boundary condition is strictly dissipative for S. This implies that the 
uniform Lopatinski condition is satisfied. 

We consider frequencies C = PC with ( close to (—1, 1, 0) where Hq(() = 
has a double eigenvalue. More precisely we consider frequencies 

(6.28) C=(-p + A,p,/> 2 7). 

In this case, we see that G is a function of a = 7 + if and p, holomorphic 
in a, as well as V, P, H and c. Moreover 

(6.29) H(() = -p 2 (aA + A~ l B + 0{p)) = p 2 H(&, p). 
The model operator is 

H(a,0) = -aA-A- 1 B=(~ a ~ 1 

\ a a + I 

H(1,0) has one eigenvalue with positive real part, with eigenvector t (b,l) 
with b = (2 + y/4 — a 2 ) /a (Note here the importance of the assumption 
a ^ 0). Therefore, for a close to 1 and p small , the negative space of 
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H(a,p) is generated by t (b(&, p),l) where b is smooth and holomorphic in 
a and b(l, 0) = b. Moreover 

(6.30) 9»6(1,0) = -(1+ 7 JL= )^o. 

a V4 - a 2 

Comparing with (|6.24|) . we see that the stability condition is violated at £ 
given by (|6.28fl , if and only if 

(6.31) b(a,p)=c(() = c(d,p). 

Given a e]0, 1[, we choose c = b and T such that the inviscid boundary 
condition reads (|6.20|) . Note that c(<7, 0) = c for all <r. Thus the equation 
Q6.31JI holds at a = 1 and p = 0. Moreover, with ()6.30|) . the implicit function 
theorem shows that for p > small, there is a(p close to 1 solution of (|6.31|) . 
providing frequencies £(p) = O(p) with 7(p) ~ p 2 > 0, where the stability 
condition is violated. □ 

4. Smooth symmetrizers. We briefly discuss here the existence of smooth 
symmetrizers for the hyperbolic operator H (|3.16|) . In the present case, we 



deduce from (|6.22j) that in polar coordinates ( = pQ, there holds 

(6.32) H((,p) = -aA + p(a 2 -fi 2 )AB + 0(p 2 ), a = ^ + i{f + f,). 

Fix £ = (1, —1,0), which corresponds to a multiple root of the hyperbolic 
part. Then a = 0, and near (£, 0) 

(6.33) H (C, p) = -A(*Id + P/3(C)£) + 0(p 2 ) 

with /?(C) = 1- Dropping the " , and changing pb to p, the matrix H is a 
perturbation for (a, p) close to (0, 0) of the following canonical example 

(6.34, (J + °), Re.>0,p>0, 

Note that (|6.29)1 derives from (|6.33j) choosing a = pa. 

Denote by E~ the negative space of H for Recr + p > 0. On can check 
directly on this example that the negative spaces have no limit as (cr, p) — > 
(0,0): the limits are different when p = and a = 0, since the positive 
spaces of A and AB are different when a^0. 

On the other hand, blowing up once more the local coordinates near 
that is taking polar coordinates (cr, p) = r(a,p), is is clear from (|6.33|) that 
E~ is a smooth function of (o",p). 
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If £(C, p) is a smooth symmetrizer for H, then (|3. 17|) implies that £ = 
X(C, 0) must be a symmetrizer for — (aA + pAB) for all a and p, equivalently 
that S = XL4 is a symmetrizer for (|6.34|) , that is 

(6.35) 5 = S* > 0, SA = AS, Re (SB) > 0. 

The first two conditions are satisfied if and only if S is diagonal and positive. 
Multiplying it by a positive factor, it must be of the form (|6.26j) . 
The third condition holds if and only if 

s > a 2 (l + s) 2 /4. 

Denoting by s m i n (a) < 1 < s max (a) < oo the roots of the equation 4s = 
a 2 (l + s 2 ), the condition reads 

(6.36) Smin(o) < s < s max (a). 

This shows that the choice of symmetrizers is much more limited in the 
viscous case compared to the inviscid one. In particular, when a is close to 
1, (|6.36|) forces to choose s in a small interval around 1. 

The boundary condition ()6.25l) is strictly dissipative for S, then (|6.2Uj) 
is strictly dissipative for S. This holds if and only if s > |c| 2 . Therefore: 

There is a smooth symmetrizer £(C, p) for H on a neighborhood of ((, 0), 
adapated to the boundary conditions (|6.25|) only if 

(6.37) |c| 2 < s max (a). 



7 The high-frequency analysis 
7.1 The main high-frequency estimate 

This section is devoted to an analysis of uniform maximal estimates for 
high-frequencies. We still assume that the Assumptions of Section 2 are 
satisfied and we prove that the anticipated (|2.58f) are satisfied when the 
uniform spectral stability conditions are satisfied, under the following ad- 
ditional structural assumptions which strengthens (H3): it means first that 
the block L 11 is hyperbolic with constant multiplicity with respect to time, 
and second that it is totally incoming our outgoing. 
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Assumption 7.1. (H5) For all u £U* and £ G M d \{0}, the eigenvalues of 
A (u, £) are real, semi-simple and have constant multiplicities . 

(H6) L n (u, d) is also hyperbolic with respect to the normal direction dxd- 

For Navier-Stokes and MHD equations and in many examples L 11 is a 
transport field 



(7.1) L n =d t + Y,a j {u)d j 

j=i 

and the condition reduces to Orf(u) 7^ for u £ U* , that is to Assumption 12.61 
which means inflow or outflow boundary conditions. The hyperbolicity con- 
dition (H6) in the normal direction is important as shown on an example 
below. On the other hand the constant multiplicity condition (H5) is more 
technical, and could be replaced by symmetry conditions: this is briefly 
discussed in Remark 17.121 

We consider the linearized equation (|2,21|) : 

(7.2) d z u = g(z,()u + f, T(()u(0)=g 

with u = \u\u 2 ,u 3 ), / = *(/\/ 2 ,/ 3 ), r as in dZSni) and g = '(g 1 ^ 2 ^ 3 ). 



Theorem 7.2. With assumptions as indicated above, assume that the uni- 
form spectral 
are p\ > a 
(f?~2l satisfy 



form spectral stability condition is satisfied for high frequencies. Then there 
are p\ > and C such that for all C, G with \(\ > pi, the solutions of 



+ +A||n 2 || i2 + ||u 3 || L 2 

+(1 + 7) V(0)| + A V(0)| + A- V 

, '" > ^cdl^ll^ + ll/^l^+A- 1 !!/ 3 !!^) 

+ C({l + i)*\g x \ + hk\g>\+ hT*\tf\). 

High frequencies require a particular analysis for two reasons. First, the 
splitting hyperbolic vs parabolic is quite different in this regime and sec- 
ond the conjugation operator <3? of Lemma 12.101 is not uniform for large £. 
The analysis is made in |MeZulj for full viscosities and Dirichlet boundary 
conditions. For partial viscosities and shocks, that is for transmission con- 
dition, the problem is solved in (JJMWZ4 . The presentation below is more 
systematic and allows for more general boundary conditions of the form 
lt2~T7T|) . 
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We now explain the general strategy of the proof. We use the notations 
\\u\\ sc = (1+7)11^11^ + A||n 2 || i2 + ||n 3 || i2 , 

ll/lls C =ll/ 1 |lL 2 + ll/ 2 ||L 2 +A- 1 ||/ 3 || i2 , 

' " ' j |«(0)| M = (1 + 7) V(0)| + A5|n 2 (0)| + A-5|n 3 (0)|, 



\g\sc = (1 +7)^l5 1 | + A2| 5 2 | + A 2\g 



3 



1) The main step in the proof of the theorem is to separate off the 
incoming and outgoing components of u. This is done using a change of 
variables u = V _1 (z, C) u which transforms the equation Q7.2JI to 

(7.5) d z u = G(z,()u + f, t(Qu(0)=g. 

There are norms similar to (|7.4|) for u and / as well; with little risk of 
confusion, we use here the same notations. An important property is that: 

(76) NU < C\\u\\ sc , ll/L < C||/||' sc , 

K0)| sc < C|tt(0)|. c , \u(0)\sc < C\u(0)\ sc , 

with C independent of (. Moreover, r(£) = r(£)V(0, C) satisfies 

(7.7) |f (C)fi(0)| flc < C|u(0)|. c . 
The new matrix Q has the important property that 

(7.8) s=(« + g °_) + g' 

with 

(7-9) \\g'u\\' sc <e(()\\u\\ sc 

where e(C) tends to as |C| tends to infinity. The block structure corresponds 
to a splitting u = (ii + ,u~) with u~ € C Nb and u + G C N+N 6 denoting 
the incoming and outgoing components respectively. 

2) One proves separate estimates for the incoming and outgoing com- 
ponents: 

(7.10) \\u+\\ sc + |fi+(0)| < C\\(8 Z - <?+)n + || sc , 

(7.11) ||u-|| sc < C\\d t - g-)u-\\ sc + C\vT(Q)\, 
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with C independent of £. (The norms are defined, identifying u G C b to 
(0, u~) G etc). As a result, with l|7.9|l . this implies that if u is a solution 
of dZSJ), then 

(7.12) \\u + \\ sc + |u+(0)| < C||/|| sc + e(0\\u\\sc, 

(7.13) \\u~\\ sc < C\\f\\ sc + e(C)IHU + C|fi-(0)|, 

3) We show that the estimates above imply that if the uniform spectral 
stability condition is satisfied, then the solutions of (|7.5|) satisfy for |£| large 
enough 

(7.14) |HU + \u(0)\ sc < C(\\f\\ sc + \g\ sc ) 
implying that the solutions of (|7.2|) satisfy 

(7.15) \\u\\ sc + |«(0)| sc < C(||/|| sc + \g\ sc ) 
that is (IT3I) . 

• Indeed, by definition, h G E~(C) if and only if there is u solution of 
d z u = Qu with u(0) = h. The corresponding u = V _1 u satisfies by (|7.13j) 

\\u-\\ sc < C\u-(0)\ +e(C)||« + ||. c 

if £ is large enough. Therefore, (|7.12j) implies that for £ large and all h £ 
E-(C), h = V _1 (0, ()h = (h + , h~) satisfies 

(7.16) \h + \ sc <e(C)\h-\ sc . 

• In addition E _ (C) := V _1 (0, f)E~(£) has dimension equal to N^, as 
the space of the h~ . Therefore, ()7.16j) shows that for £ large, the projection 
h i — ^ h is bijective from E~(£) to C Nb , with inverse uniformly bounded in 
the norm | • | sc . 

The uniform spectral stability condition reads 

(7.17) V/i G E-(C), \h\ sc < c\r(Qh\ sc 
(see ()2.59j0 . Using (|7.6|) . this implies 

(7.18) vkr(C), \h\ S c < c\r(c)h\sc 

Using the isomorphism between E~(£) and C Nb , we see that for £ large 
enough and h~~ G C Nt >, there is h + such that (h + ,h~) G E~(C). Together 
with (f?~T!)|) and ({777jl . there holds 

< I^U < C\f(()h\ sc < C\t (C)(0,fe-)| ac + e(C)|fe _ |- c . 
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For £ large, the last term can be dropped, increasing C. Finally, we conclude 
that for all h £ C N 

(7.19) \h\sc<C\f(()h\ sc + C\h+\ sc . 

Applying this estimate to u(0), combining with (|7.1U|) and (|7.11|) and ab- 
sorbing the error term Q'u for £ large, we immediately obtain (|7.14|) . 

The third part of the proof will not be repeated. We will focus on the 
reduction (|7,5|) and on the proof of the estimates for u . 

7.2 Spectral analysis of the symbol 

Consider the linearized operator (|2.2U|) 

-Bd 2 + Ad z +M. 

The coefficients satisfy 

B{z) = B dd {w(z)) 

d-i 

A{z, C) = A d {w(z)) - iVj {B jd + B dJ ) (w(z)) + E d {z) 

3=1 

d-1 

M(z,C) = (ir + j)A (w(z)) +^2ir h (A j (w(z)) + E 3 {z)) 

i=i 

d-1 

+ XI VjVkBj, k {w{z)) + E (z) 
j,k=i 

where the E\~ are functions, independent of C, which involve derivatives of 
w and thus converge to at an exponential rate when z tends to infinity. 
Moreover, we note that 

(7.21) El 1 = 0, El 2 = for k > 0. 

With (|2.2j) . we also remark that M 12 does not depend on r and 7 

We start with a spectral analysis of the matrix Q in (|2.21j) . It is conve- 
nient to use here the notations u = (u 1 ,u 2 ,u 3 ) E (£ N ~ N x x . In 
the corresponding block decomposition of matrices and using the notations 
above, there holds 



(7.20) 





(Q 11 


Q 12 


g 13 


(7.22) 


g= 





Id 




W 1 


g32 


g33 
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where 



G n = 


-(A n y 




g31 


= (B 22 ) 


-\A 21 g n 


+ M 21 ) 


gu = 


-(A n y 




g32 


= (B 22 ) 


-\A 21 g 12 


+ M 22 ) 


g 13 = 


-(a 11 : 


rU 12 , 


g33 


= (B 22 ) 


~l {A 21gV3 


+ A 22 ). 



Note that g , G 12 -, g 31 and g 33 are first order (linear or affine in £), that 
g 32 is second order (at most quadratic in £) and that g 13 is of order zero 
(independent of £). We denote by g® b their principal part (leading order 
part as polynomials). We note that 

(7.23) g« b (z,() = Gf(w(z),() when (a, b) ^(3,1), 

with 

GlHuX) = -(AYiu))- 1 ^ + ir)A^(u) + J2 d J - 1 li 7 lj AY(u)), 

Gl 2 (uX) = -{Af{u))-^U^ A f^ 
Gf{u) = -{A\\u))^Af{u) 

Gf{uX) = (B^y'EtkliViVkB 22 ^)), 

Gf(u, C) = -(B 22 (u))~ l Y,U *<* {B 2 %u) + Bg(«)) . 

The principal term of g 3,1 involves derivatives of the profile w. Denoting by 
p = lim^+oo w(z) = w(oo) the end state of the profile w, we note that the 
end state of g 31 is 

SJV.C) = {B 22 {p))- l {{ 1 + ir)A 2 \p) + Y,ir ]j Af (p) + Af(p)G?(p,0). 

3=1 

There are similar formulas using the matrices Aj and Bj^ of (|2.3j) . 

The spectral analysis is easier when all the terms are reduced to first 
order. If u = (u l ,u 2 ,u 3 ) is replaced by u = hi^i-u := {u l ,u 2 ,\C,\~ l u 3 ), g is 
replaced by 




(7.24) g = h lcl gh^=[ KIM ] := ( £1 £22 



with obvious definitions of V ab . Note that g is or order one, while V 21 is of 
order zero. Thus 

(7.25) g(z, o = g P (z, o + o(i), <? p = ^ ^2) = o(ICI). 
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Moreover, since the coefficients in Q converge exponentially at infinity, the 
remainder in (|7.25|) is uniform in z G M+ and \(\ > 1. Moreover, the principal 
part of V 22 is if the form Vf{z, Q = P 22 (w(z), (). 

Lemma 7.3. i) For all ( G with 7 > and f] ^ and for all and 

z > 0, Gr>{z, C) has no eigenvalues on the imaginary axis; moreover, the 
number of eigenvalues in {Re n < 0} is = N\ + N' . 

ii) for all compact subset ofU*, there are c > and 5 > such that for 
all u in the given compact and all ( G such that either 7 < 5\(\ or 

\t]\ < o~\C\> the distance between the spectrum of G^(u, £) and f/ie spectrum 
of P 22 (u,() is larger than c\(\. 

Proof. The spectrum of Q p is the union of the spectra of Gp 1 and P 22 . By 

homogeneity, it suffices to consider ( G S^. 

a) Gp 1 is related to L 11 since A^(i^ + G^-(u, Q) = L u (u, 7 + ir, irj, 
By Assumption (H3), L 11 is hyperbolic in the time direction, hence Gp 1 has 
no eigenvalues on the imaginary axis when 7 > 0; moreover, the boundary 
is noncharacteristic for L 11 by Assumption 12.61 implying that the number 
of eigenvalues of Gp 1 in {Re fi < 0} is equal to the number of positive 
eigenvalues of A\ , that is is iVi. 

Next, note that 

p22 _ ( |C|Id\ 

^ -\\C\~ l Gf Gf) 

Thus, i£ is an eigenvalue of P 22 if and only if is an eigenvalue of B 22 (r], £), 
which is impossible by (H2) if 77 7^ 0. Thus, the eigenvalues of P 22 are 
not purely imaginary when 77 7^ 0. Moreover, the number of eigenvalues in 
{Re/i < 0} is N' (see |MeZulj ). This finishes the proof of i). 

b) If 7] = 0, Gp 2 and Gp 3 vanish, hence the spectrum of P 22 is {0}. On 
the other hand is not an eigenvalue of Gp 1 = — (7 + ir)(A^-) Aq since 
Aj- and Aq 1 are invertible and I7 + ir| = \C\ = 1. 

If 7 = and r\ 7^ 0, the eigenvalues of P 22 are not in iM.. On the other 
hand, by Assumption (H6) the eigenvalues of Gp 1 are purely imaginary, thus 
P 22 and Gp 1 have no common eigenvalue. This finishes the proof of ii). □ 

The analysis in a purely "elliptic" zone {7 > <5|£| and \r]\ > d\C\} with 
5 > 0, is easy, see below. The most difficult and important part is to 
understand the "hyperbolic-parabolic" decoupling in an arbitrarily small 
cone 

(7.26) C 5 = {0< 7 <<5|C|}U{|r ? |<5|C|} 
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with 5 such that property ii) of Lemma [731 holds for u in a simply connected 
neighborhood Uq of a compact set which contains the curve {w(z), z E 
[0, +oo[}. There, the usual homogeneity and the parabolic homogeneity are 
in competition, leading to different classes of symbols. We use the following 
terminology: let £ = (t, 7,7/) and for a multi-index a = (a T , 0^,07) E 
N x N d_1 x N, set 

|a| = a T + |a,j| and (a) = 2(a T + a 7 ) + |a^|. 

Recall that the parabolic weight is A = (1 + r 2 + 7 2 + 1 77 1 4 ) 4 . 

Definition 7.4. i ) r m (0) denotes the space of homogeneous symbols of 
order m, that is of functions h(z, Q E C°°(M+ x ft) such that there is 8 > 
such that for all a E N rf+1 and a// k E N, i/iere are constants C a ^ such that 
for Id > 1 : 

(7.27) \d?h\ < C a ,o|Cr- |a| , i/* = 0, 

(7.28) < C Qifc e- fe |Cr- |cl1 , i/ fc > 0, 



ii j Pr m (r2) denotes the space of parabolic symbols of order m, that is 
of functions h(z,() E C°°(1R+ x f2) satisfying similar estimates with \( i \ m ~\ a \ 
replaced by A m ~^ . 

VKe use t/ie same notation for spaces of homogeneous or parabolic matrix 
symbols of any fixed dimension. 

Lemma 7.5. For all ( E S d Pi C5, i/iere is a a conical neighborhood U o/( 
and t/iere are matrices Wp 2 E r°(Jl) and Wp 1 , homogenous of degree in £ 
/or u E ZYq suc/i £/iai 

(7.29) W 21 ^ 1 - P 22 W 21 = ICIPj 1 . 

(7.30) ^p 1 ^ 2 - YJfvf = -^ 2 . 

Proof. By homogeneity, it is sufficient to construct Wp 31 for |£| = 1. By 
Lemma E21 for ( E S d+1 n Ca and u e W *, the spectra of G^(u,C) and 
P 22 (u, £) do not intersect, so that the linear system of equation 

XG^{u,Q-Pl 2 {u,C)X = Y 

has a unique solution X = X(u,£)Y. Therefore W 21 (z, () = \(\ X(w(z), ()V 21 (z, () 
satisfies RUM (Note that V 21 is of degree 0). 

The construction of Wp 2 is similar, noticing that "Pp 2 is of degree 1. □ 
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In the block structure of G, there holds 

(7-31) <=(lf)' W p 2 = ( V p 2 V p 3 ) 

and (f?~2T)|) reads 

(7.32) Vp^-KlVf = 0, 



(7.33) Vfgf -\Q\~ l Gfvl l -Gfvf = Gf 



Similarly, 



(7.34) G^x 2 -\crxx = -% 



(7.35) £ii v i3_| C | V i2_vi 3 Sf = | C |gi3. 



v 12 


icr^ 3 


Id 








Id 



For further use, we make the following remark : by (|7.23|) . we see that Gp 2 
and Gp 2 vanish when r\ = 0. Therefore, ()7.34|) implies that V 12 also vanishes 
when T) = and hence 

(7-36) V 12 (z,C) = 0(N/|C|). 

With these notations, let 

/Id 0' 

\ V 31 Id/ \0 

and V = V/V//. Using the conjugation u = Vu, f = Vf, for £ in the in the 
cone Cs, the equation ()7.2[) is transformed to 

(7.37) d z u = Gu + f, tu(0) = g 

with G = V^GV - V~ l d z V and f (C) = T(C)V(0, C). 
Lemma 7.6. The entries of G satisfy: 

G 11 - {G u + \cr l G 12 vl l + s 13 v 31 ) e r- 1 , 

£12 e pO^ £13 g p i £21 g £31 g pO^ 

G 22 er°, g 23 -Ugt-\ 

£32 _ (£32 _ y31gl2) g pO^ £33 _ g33 £ pO^ 
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Proof. We first compute the entries of Qi = Vj 1 QVi. Direct computations 
show that 

gf = gii + |crV 2 v 21 + e 13 v 31 , gf = g 12 , Qf = g 13 

Moreover, 

Q f = -icrx's 11 + v 31 - icrv^icrwf +g i3 v p 31 ). 



The first two terms are of degree zero, and by (|7.32j) . the sum of their 
principal terms vanishes; the third term is of degree —1 thus Qf E T -1 . 
Similarly, Qf is of degree 1 and its principal part vanishes by 1)7.33(1 . Thus, 

Gf e r~\ Qf g r°. 

Next 

gf = _| C |-i Vp 2 V 2 g r°, Qf - id = -icr'v 21 *? 13 g r- 1 . 

The computations for Qu = VJt QiVn are quite similar. This new con- 
jugation annihilates the principal parts of Qf and Qf and contributes to 
remainder terms in the other entries. 

Finally, direct computations show that V d z V only contributes to re- 
mainder. □ 

The main idea is to consider (|7.37|) as a perturbation of the decoupled 
system 

(7.38) d z u l = Q 11 u 1 + h, 

*(?)-(££)(?)♦(? 

Introduce then 





(Q 11 








(7.40) 


g' = g-\ o 





Id 




V o 


£32 


£33 



The next lemma how the estimates are transported by the change of variables 
u = Vu. We use the notations (|7.4|) for the scaled norms. 
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Lemma 7.7. There are constant C and p\ such that for all £ in the cone 
C$ with \(\ > p\, there holds 



\V- l u\\ sc < C\\u\\ sc , \\Vf\\' sc < C\\f\\' sc , 

V^mUc < C\u(0)\ sc , |Vn(0)| sc < C\u 



I sc- 



(7 ' 41 > IV-- 

and 

(7.42) |f (C)«(0)| sc < C\u( 
Moreover, 

(7.43) \\Q'u\\sc < CA _1 ||u| 
Proof. Direct computations, using (|7.36f) . show that u = Vu satisfies 

u 1 = o^u 1 + od^i icrv + oacrv, 

U 2 = 0(|Cr 1 )ti 1 + 0(l)^ + 0(|C|- 1 )n 3 , 
u 3 = O(l)^ 1 + 0{l)u 2 + 0(l)u 3 . 

This implies the first estimate in (j7.41[) . using the inequalities 

(1 + 7)H/ICI<A, (l + 7 )/|C|<l, A/|C| < 1- 

The proof of the other estimates of (|7.41|) is similar, using in particular for 

the traces the inequality (1 +7) 5 |?7|/|C| ^ A.2. 

The inequality (|7.42|) follows from the second line of (|7.41|) and the 

estimate |ru(0)| sc < | 

^(0)|sc which is a direct consequence of the form (|2.56j) 
of the boundary conditions. 

Finally, Lemma 17.61 implies that / = Q'u satisfies 

; 2 = o(icr v + o(i)£ 2 + oacrv, 

f = 0(l)v} + 0{l)u 2 + 0(l)u 3 . 
and (TTIal follows. □ 



The parabolic bloc (|7.39f) is studied in [MeZulj . We now focus on the 
hyperbolic block (|7,38|) . recalling and extending the analysis of |(tMWZ4] . 
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7.3 Analysis of the hyperbolic block. 
7.3.1 The genuine coupling condition 

For u €U*, denote by Xj(u, £) the distinct eigenvalues of A U (u, £), which 
are real and have constant multiplicity Vj by Assumption (H5). Assumption 
(H6) implies the following: 

Lemma 7.8. For all u £ U* , all £ 6 R d and a// j, i/iere holds <9,e d Aj(u, £) / 
0, and all these derivatives have the same sign. 

Proof. If d^ d Xj(u, i], £, d ) = 0, then the equation r + A(n,£d) = would have 
complex roots in ^ for some r close to r = —Xj(u,r],^ d ) (recall that Xj is 
real analytic) . Thus hyperbolicity in the normal direction prevents glancing. 
Moroever, by continuity the sign of d^ d Xj(u, n, £d) is constant for all 6 R 
when 77 7^ 0. Thus the functions ^ i-> Aj(u, n,£d) are monotone and tend 
to infinity as ^ tends to ±00. Since Xj / A& when j k, they must be all 
increasing or all decreasing. This remains true for n = by continuity. □ 

According to the terminology of Section 4, we will say that the hyperbolic 
block L n is incoming [resp. outgoing] when the derivatives d^ d Xj(u,£,) are 
positive [resp. negative]. 

Corollary 7.9. i) The matrix GV-(u,C) has no purely imaginary eigenvalues 
when 7 > 0. They are all lying in {Re fi > 0} if the 11-block is outgoing and 
in {Re fj, < 0} if it is incoming. 

ii) Near points £ with 7 = 0, Gp 1 (u, C) semi-simple eigenvalues 
/ij(n, C) 0/ constant multiplicity Vj, which are purely imaginary when 7 = 0. 
Moreover, d^Refij > when the 11-block is outgoing and d^Refij < when 
the 11-block is incoming. 

Proof. Note that fj, is an eigenvalue of Gp^u, C) if and only if — r + 27 is an 

eigenvalue of A U (u, n,£) with £ = — i/i. 

Consider the equations in ^ : r + Xj(u, n, £d) = 0. Since A-,- is strictly 
monotone and tends to infinity at both infinity, it always have a unique 
solution, ifjj(u,r],T) and Qj-V'j has the same sign as —d^ d Xj. This solu- 
tion extends analytically for Imr small. This yields distinct eigenvalues 
Hj(u,Q = iipj(u,rj,T — ij) of Gp 1 for ( close to the real domain. In partic- 
ular d^fij = d T tpj and the eigenvalues all lie in {Refi > 0} if the 11-block is 
outgoing and in {Re jjl > 0} if it is incoming. 

The kernel of Gp 1 — \ij is the kernel of A — Xj , thus has dimension equal 
to the multiplicity of Xj. Since these dimensions add up to N 1 , this shows 
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that Gp has only semi-simple eigenvalues of constant multiplicity, which all 
lie in a given half space when 7 > 0. 

Hyperbolicity of L 11 implies that Gp (u, has no purely imaginary 
eigenvalues when 7 / and by continuity they all lie in the same half 
space. □ 

Next we need more information on the zero-th order correction of Q . 
From GHM (17311 and (t?~22]) we see that 

(7.44) g 11 (z,()-(V~ 1 d z V) 11 )=G 1 p 1 (w(z),()+£(z,0, 
where £ £ T°. Denote its principal part by £ p . Its limit at z = 00 is 

(7.45) E p (p,0 = \C\- 1 g 1 p 2 (p,c)v p 21 (p,c)g1 3 (p,C)v p 31 (pX) 

where p = lim z ^ +00 w(z) and V p 21 (p, £), V p 31 (p, () denote the end points of 
Vp 1 and Vp , that is the solutions of the intertwining relations (|7.32|) (|7.33|) 
with matrices G p b replaced bay their endpoint values G p b (p, (). The next 
result is crucial and follows from the genuine coupling condition (H4). 

Proposition 7.10. Fix C with \(\ = 1 and 7 = 0. For £ in a neighborhood of 
consider a basis where G n (u, £) has the block diagonal form diag(//jld^ ) . 
Denote by Ej^{u,C]) the corresponding blocks of E is this basis. Then, for 
u £ U the eingenvalues of the diagonal blocks Re Ejj have a positive [resp. 
negative] real part if the 11-block is outgoing [resp. incoming]. 

Proof. It is sufficient to prove the positivity at £. Suppose that 7 = 0, denote 
by ipj, p with pG {1, . . . , Uj} a basis of eigenvectors of G n (u, £). Fix j and 
set £d = —i/j,j(u,C) £R, f= (Vi^d)- Then the cpj tP are right eigenvectors of 
A (it,£) associated to the eigenvalue — r = Xj(u, £). 

Consider left eigenvectors £j )P of A (it, £), dual to the ipj tP . Then, the left 
eigenvectors of Gp 1 (u, () associated to fij are -^-ljA d with f3j = d^ d Xj(u, 77, £, 
wee Lemma 14.191 The entries of the block Ejj are 

(7.46) —i^AfEpiuX)^. 

Pj 

Computing the eigenvalues of order e of B(u,£) + ieA(u, £), leads to consider 
the matrix 

(7.47) iA U +eA 12 (B 22 )- l A 21 . 
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The genuine coupling condition (H4) implies that for u 6 U, its spectrum 
lies in Re p > ce for e small, and this implies that the matrix Fjj with 
entries 

(7.48) t j;P A 12 (B 22 )- 1 A 21 <p j;P , 

has its eigenvalues in the right half plane {Re/i > 0}. 
Because Gp<pj ;P i = i£d ( Pj,p', the relation (|7.32f) implies 

= \c\-X 1g p^p' = ^icrv p 2 v iy , 

and, using the expressions of the matrices G a,b yields 

(icr'Gix'+Gix 1 )^' = -i\c\- 1 <& 1 r r z a (.v,€)vZ L <Pj*> 

and 

By (|7.33f) this is equal to 

-GpVjy = -^f^A 21 ^^)^. 

Thus 
and 

Multiplying on the left by IjA^, this shows that the coefficients in (|7.46JI 
and ()7.48[) only differ by the factor —l/(3j, and the proposition follows. □ 

7.3.2 Estimates 

We are now in position to prove maximal estimates for the solutions of the 
equation (|7.38|) . 

Proposition 7.11. There are constants C and p\ > 1 such that for all £ 
in the cone C$ with \(\ > p\ and all u 1 and f 1 in L 2 (R + ) satisfying l|7.38jl . 
there holds 

(1 + 7 )||« 1 ||l2 + (1+ 7 )^ 1+ (0)| 

( 7 - 49 ) / -1 1 1 , \ 

<C7(||/ 1 || ia + (l+ 7 )3|«i-(0)|) 

where -u 1+ = u 1 and u l ~ = if the ll-block is outgoing and u l+ = and 
u 1 ^ = u 1 if it is incoming. 
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Proof, a) Fix £ E S + . We prove the estimate for £ in a conical neigh- 
borhood of C- Suppose first that 7 = (the most difficult case). By 
Corollary 17.91 there is a matrix V 11 (z, () homogeneous of degree such that 
(V 11 ) -1 £/p 1 V 11 = diag(/j,j(w(z), C)Id Vj -. Setting u 1 = V 11 ^ 1 transforms the 
equation to 

(7.50) dzu 1 = (diagO^Xz), C)Id Vj . + Sju 1 + / x 

with £ = £ — (V 11 )~ 1 d z V 11 € r , whose principal part £ p has the same end 
point Ep(p,C) as £ p . 

As usual, since the fXj are pairwise distinct, there is a new change u 1 = 
(Id + V-i)^ 1 with Vl\ G r _1 , such that the resulting system has the same 
form with the additional property that the zero-th order part is also block 
diagonal, so that £ p = di&g(£jj) and the end points of the blocks £jj are 
Ejj introduced in Proposition 17.101 

The term {£ — £ p )u is 0(|C| _1 |^|), is incorporated to J 1 and finally 
absorbed from the right to the left of the inequality by choosing |£| large 
enough. This reduces the proof to the case where the equation reads 

(7.51) dzu 1 = fij(w(z), On 1 + ^(C^ 1 + Fjjiz, On 1 + f 

with \F jtj \ < C e- 9z . 

Consider the outgoing case. Then, Corollary 17.91 implies that there is 
a constant c > such that He [Aj(u,Q > cry. Moreover, Proposition 17.101 
implies that the eigenvalues of Ej have a positive real part. Thus, there is 
a positive definite (constant) matrix geld such that Re SEjj is definite 
positive, say ReSEjj > Id. Introduce a = Cq\S\ Jq e ds such that d z a > 
\SFjA and a is bounded in L°° uniformly with respect to (\ Therefore, 
multiplying the equation by e 2a ^ S and taking the L? scalar product with 
u l implies that 

(1 + c 7 )||eV||i 2 + |^(0)| 2 < C\\e a u 1 \\ L2 \\e a f 1 \\ L 2 

which implies (|7.49|) . The proof in the incoming case is similar. 

b) Suppose next that 7 = 0. Consider again the outgoing case. Then, 
the eigenvalues of Gp 1 satisfy Re/i > c|£| in a conical neighborhood of £. 
This is the classical "elliptic" case. There is a symmetric definite positive 
matrix S(u,() S T such that ReSG 11 > c|C|Id and usual integrations by 
parts imply that 

clClll^llis + |^(o)| 2 < ciiu 1 !!^!!/ 1 !!^ + c^wl* 
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where C\ involve estimates of the zero-th order terms, which include d z S(w(z), (). 
This term is eliminated choosing \Q large enough. The proof in the incoming 
case is similar. □ 

Remark 7.12. The proof above contains two ingredients. First, the 11- 
block is totally incoming or totally outgoing, in analogy with the terminol- 
ogy of Section 4. Thus the decoupling incoming/outgoing is trivial. More 
generally, this could be replaced by a decoupling condition in the spirit 
of Section 4. For instance, for shocks, such a decoupling is immediate in 
GMWZ4 corresponding to equations on each side of the front. Next, we 
construct symmetrizers for the incoming and outgoing components. There 
we use the genuine coupling condition. If the eigenvalues are not of constant 
multiplicity one can introduce adapted bases or use symmetry also in the 
spirit of Section 4. 



7.3.3 About Assumption (H6) 

We show on an example that hyperbolicity in the normal direction is crucial 
in the proof of estimates of the form ()7.49j) . Suppose that the L 11 - block 
reads 



(7.52) 



dtu — d y u + d x v, 
dtv + d y v + d x u. 



Then, on the Fourier side, the 11 equation will be of the form 
(7.53) 



(i(r - 77) + j)u + d z v + a(z)u = f, 
(i(r + r]) + j)v + d x u + a{z)v = g, 

and the only information we have from the genuine coupling condition is 
that a is positive at z = +00. Suppose that a(z$) < for some zq > 0. 
Then glancing waves for (|7.52j) will propagate parallel to the boundary and 
thus may remain in a region where a is negative and thus may never be 
damped. This is illustrated by choosing r = 77, large, 7 = — a(zo) > and 

u t(z) = x{t~H z - z o)), v T (z) 



2ir + 7 + a 



with x G Co° W- Then (T7331 is satisfied with / = (a(z) - a(z ))u T + d z v T 
and g = 0. Moreover, [|/[|i2 = 0(r - 3)||it||£,2 and u(0) = v(0) = 0, showing 
that no estimate of the form (|7.49|) can be valid. 
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7.4 Proof of Theorem \77Z\ 



7.4.1 In the cone C5 

We consider now the equation (|7.39j) and briefly recall the results from 
[MeZulj . It is natural to rescale the problem using the parabolic weights: 
with v 2 = u 2 and u 3 = A -1 -u 3 and g 2 = f 2 and g 3 = A -1 / 3 the system reads 



(7.54) 94$=g P K) + K 



with 



of quasi- homogenenous degree one and principal part Gp(w(z),Q with 

(7 - 55) GPM) = (a-H^ + t)^)- + C? M ) Gfltv)) 

Lemma 7.13 ([MeZulJ). There is c > such that the spectrum of Gp 
lies in {|Re//| > cA}, with N' eigenvalues, counted with their multiplicity, 
of positive real part. There is a smooth change of variables W G Pr° such 
that 

with V± G Pr 1 having their eigenvalues satisfying ±Re// > cA. 
Introduce 

^ - w- 1 (i 



V J \v 

Corollary 7.14 ( [MeZulJ). There are C and p\ such that for all £ £ 
with I C| > pi, there holds 

A\\v + \\ L 2 + A3| v +(0)| < C\\(d z - V + )v + \\ L 2, 
< C\\(d z -V~)v-\\ L 2 +C7A5|«-(0)|. 
Scaling back, introduce 

™ (£)-(? aXo> (£>(o a)-C- 

Because, W~ 1 d z W is uniformly bounded, the Corollary implies the following 
estimate: 
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Proposition 7.15. There are C and p\ such that for all ( G Cj with |£| > 
px, there holds 

A||u 2 ' + || i2 + ||n 3 ' + || L 2 + A^|n 2 ' + (0)| + A-^|n 3 -+(0)| 

< c\\f 2 \\ L 2 + ca- 1 !!/ 3 !!^ + \\u 2 \\ L2 + CA- 1 !^ 3 !!^, 



A||u 2, "||La + ||« 3 '"||La < Chh\u 2 ^ (0)| + CA~3 |w 3 '~(0)| 

+ C\\f 2 \\ L 2 + CK^WpWv + ||^ 2 || L2 + Cfi 
Finally, with u 1 '^ as in Proposition 17. Ill introduce 



(7.57) u ± = t (u 1 ' ± ,u 2 > ± ,u 2 > ± ) 

Adding up the various estimates and using H7.43j) , one obtains the following 
estimates. 

Proposition 7.16. There are C and p\ such that for all (, £ C$ with |£| > p\ 
and all u G H 1 (W + ): 

(7.58) \\u+\\ sc + \u+(0)\ < C\\(d z - G)u\\ sc + A~ 1 ||^|U C , 

(7.59) \\u-\\ sc < C\\(d z -g)u|| sc + A- 1 ||u|| sc + C|u-(0)|. 

As indicated at the end of Section 7.1, these estimates imply the max- 
imal estimates of Theorem 17.21 provided that the boundary conditions are 
uniformly spectral stable. 

7.4.2 Analysis in the central zone 

We now consider the remaining cone where 

(7.60) ( G 7>£|Cl and |r?| > 6\(\. 

We consider the rescaled Q matrix (|7.25|) . for the rescaled unknows u = 
h, f ,« := (uSuMCI- 1 ^), / = h K |/ := (/ 1 ,/ 2 ,|Cr 1 / 3 )- We note that in 
the region under consideration we now have (1 + 7) ~ A ~ |£|, so that the 
rescaled norms (|7.4|) are equivalent to 

\\u\\sc ~ |CllNU 2 i 

(7.61) K0)| sc «|C|^R0)|, 
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By Lemma 17.31 there is a smooth matrix V £ T° such that 





tt + has dimension N + N' — Nj, and ti has dimension iVj,. The equation for 
u reads 



with = C/ dia s + 0(1)- The ellipticity of C/ diag immediately implies the 
following estimates. 

Proposition 7.17. There are constants C and p\ such that for all £ satis- 
fying (|7.6U|) and \Q\ > p\ and all u 6 // 1 (R + ) satisfying (|7,63|) . i/iere ZioWs 



Thanks to (j7,61|) . this is the exact analogue of Proposition 17.161 and 
these estimates imply the maximal estimates of Theorem 17.21 provided that 
the boundary conditions are uniformly spectral stable, as explained in Sec- 
tion 7.1. 

8 Application to magnetohydrodynamics 

We now apply our results to the equations of isentropic magnetohydrody- 
namics (MHD), for which the inviscid case was treated in |MeZu2j . The 
full (nonisentropic) inviscid equations have been treated in |Kwj . and have 
essentially the same symbolic structure as the isentropic inviscid equations. 

8.1 The equations 

The equations of isentropic magnetohydrodynamics (MHD) appear in basic 
form as 



(7.63) 



d z u = Qu + f 



(7.64) 
(7.65) 



\C\\\u + \\ L 2 + \C\*\u + (o)\ <c\\f\y + c\\u\\ L 2, 

ICI Wu'Wv < C\\f\\ L 2 + C\\u\\ L 2 + C|C|^|n-(0)| 2 . 



(8.1) 



' d t p + dw(pu) = 

< dt(pu) + div(pu*it) + Vp + H x curlff = evAu 
k d t H + curl(H x u) = epAH 



S3 



(8.2) divH = 0, 

where p E R represents density, m£K 3 fluid velocity, p = p{p) G R pressure, 
and i7 G R 3 magnetic field. When H = 0, IjH.lJ) reduces to the equations of 
isentropic fluid dynamics. We assume that v and p are positive. 
Equations (|8.1|) may be put in conservative form using identity 

(8.3) H x curltf = (l/2)div(|#| 2 / - 2#'if) tr + HdivH 
together with constraint Q8.2JI to express the second equation as 

(8.4) d t (pu) + div(pu*u) + Vp + (l/2)div(|#| 2 / - 2F*F) tr = evAu. 

They may be put in symmetrizable (but no longer conservative) form by a 
further change, using identity 

(8.5) curl(F xu) = (divu)F + (u ■ V)H - (divii> - (H ■ V)u 
together with constraint ()8.2|) to express the third equation as 

(8.6) d t H + (dwu)H + (u • V)H - (H ■ V)u = peAH. 

Forgetting the constraint equation, we get a 7 x 7 symmetric system. 

Neglecting zero-th order terms, the linearized equations of (|8,1|) about 
(p, u, H) are 



1.7) 



' D t p + pdivu 

pD t u + c 2 Vp + H x curLff - evAii 
k D t H + (divu)H - H-Vii- epAH 



with Dt = dt+u-\/ and c 2 = dp/ dp which we assume to be positive. This sys- 
tem is hyperbolic symmetric, with symmetrizer S = block-diag(c 2 , pld, Id). 
It enters the general framework of linearized equations studied in this paper, 
with parameters (p,u,H). 

8.2 Eigenvalues and eigenvectors 

Eigenvalues and eigenvectors of the symbol solve 
fa + p(£ • u) = 

fii + c 2 &^ + v x (£ x t>) = ii/|£ 2 |£t/p 
. tu + (£ • — (f • £,)u = ip\t,\ 2 i). 



84 



with 



(8.9) f = T + u-t, v = H^fp, & = p/p, v = H/y/p. 

The structure condition (|2,2|) is satisfied with N' = 6. The kernel of -B(£) is 
generated by (1, 0, ... 0) which is never an eigenvector of A(£) when £ ^ 0. 
Thus the Assumptions (HI'), (HI) and (H2) are satisfied. 

Consider next the inviscisd problem. The seven eigenvalues of A(p, u, H, £) 
arc (see e.g. |MeZu2| 1: 



U0) 



A = u • £, 
X± s = A ± c s |£|, 
A ±2 = A ± v • £, 
I A±/ = \o±c f \£\, 



c 2 f := ]-(c 2 + \v\ 2 + \J{c 2 - \v\ 2 ) 2 + 46 2 c 2 



with 

r 2 : = 

2 

i (c 2 + |u| 2 - ^(c 2 -\v\ 2 ) 2 + 4b 2 c 2 ^ 

c 2 = P '(p)>o, v = h/^p, b = \ix V \, | = e/iei. 

The first eigenvalue corresponds to the transport of the constraint. It can 
be decoupled from the system : there is a smooth one dimensional subspace, 
Eo such that A(£) = Ao on this space and Ejj- is stable for A(£). The other 
eigenvalues are in general simple. 

Lemma 8.1 ([MeZu2]). Assume thatO < \v\ 2 ^ c 2 . Consider £ £ R 3 \{0}. 
When £ ■ v ^ and ^ x H ^ 0, the eigenvalues are simple, 
ii) On the manifold ^ x v = 0, Xo is simple. When \v\ 2 < c 2 [resp. 
\v\ 2 > c 2 ], X±f [resp. X± s ] are simple, the other eigenvalues A±2 = A± s 
[resp. X±2 = X±f] are double, algebraically regular but not geometrically 
regular. Moreover, 

(8.11) X ±2 -X ±s = 0{\Zxv\ 2 ) [resp. X ±2 -X ±f = 0(\^xv\ 2 )}. 

Hi) On the manifold £ ■ v = the eigenvalues X±f are simple and the 
multiple eigenvalue Ao = A± s = A±2 is geometrically regular. More precisely, 
there are smooth X±i such that {A S ,A_ S } = {Ai,A_i}. Moreover, 

(8.12) X ±1 = u.a±5vZ + O{(v0 2 ), 6= ' • 

yr + h z 
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One can choose smooth eigenvectors eo,e±i,e±2 such that, on the manifold 
£-v = 0, 

1 / ° 

e ±l = 7=t r I Tv/d I , e ±2 = -7=— I > 





V2\v\ ' V2\ 



v 



8.3 Glancing and viscous coupling 

The boundary {23 = 0} is noncharacteristic for the hyperbolic part if and 
only if 

(8.14) u 3 i {0, v 3 , ±c.(n), ±c f (n)} 

where c s (n) and c/(s) are the slow and fast speed computed in the normal 
direction n = (0, 0, 1). 

Lemma 8.2. Assume that < \v\ 7^ c. 

i) On the manifold ^ x v = 0, the multiple eigenvalues are nonglancing 
if and only if u 3 7^ ±^3. In this case, they are totally nonglancing. 

ii) On the manifold £-v = 0, the multiple eigenvalues are non glancing if 
and only if u 3 7^ ; u 3 7^ ±v 3 and u 3 7^ ±<5t>3. They are totally nonglancing 
when 1 1*3 1 > 1 173 1 . 

Proof. By (|8.11j) . on £ x t> = 0, with j = s when \v\ < c and j = f when 
\v\ > c, there holds 

%A ±j = %A ±2 = «3 ± v 3 . 

This implies i). 

In addition, c^Ao = 7x3, <%,A±2 = ^3±^3> and by (|8.12l) <9g 3 A±i = 1*3 ±6^3 
on the manifold t;v = 0. This implies ii). □ 

Next we study the viscous coupling of vectors ej at geometrically reg- 
ular modes. In the variables (/?/ p,ii,v), the system (|8.7|) is symmetric, 
with symmetrizer S = diag(c 2 , Id, Id), and the viscosity matrix is -B(£) = 
|£| 2 diag(0, vld/p, pld). The basis (|8.13j) is orthonormal for S. Therefore, 
according to the general rule Q6.9f) . the matrix is symmetric with nondi- 
agonal entries 

-°0,±1 — n 0,±2 _ -°±1,±2 — U ' 

2 2p' 2 ~ 2 2 2p" 
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When |ti3 1 < 1 173 1 , then one of the eigenvalue A±2 is incoming and the other 
one outgoing (depending on the sign of v 3 ). Therefore, if p, — v/p, the 
coupling coefficient B\ _ 2 does not vanish. Summing up, we have proved: 

Lemma 8.3. If \uz\ < |f 3 1 , and v 7^ pp, then the decoupling condition (|4.9[) 
is not satisfied at modes where £ ■ v = 0. 

Remark 8.4. The decoupling of the mode Ao from the other ones reflects 
that the constraint (|8.2j) is propagated by the viscous equation as well. The 
other partial decoupling observed above depend on the particular choice of 
the viscosity matrices and disappear for general B. 

8.4 Shocks 

Consider an inviscid planar shock. We suppose that the front is x 3 = at 
and denote by (p~~ ,u~ , H~) and (p + ,u + ,H + ) the states on the left and on 
the right respectively. All the analysis of the preceding section is valid, if 
we change u 3 to u 3 — a. 

The jump conditions are deduced from the conservative form of the equa- 
tions: 



where r 3 = *(0,0, 1). The last jump condition comes from the constraint 
equation ()8.2j) . Apparently this system of 8 scalar equations is too large. 
However, projecting the third equation in the normal direction yields o-[H^\ = 
which is implied by the last equation. This shows that (|8.16|) is made of 
7 independent equations, as expected. Denoting by n tg and H tg the tan- 
gential part of u and H, that is their orthogonal projection on (|8.16j) is 
equivalent to 



(8.16) 



< 




(8.17) 



' [P(U3 ~ -7)] = 0, 

< [pu(u 3 - a)} + r 3 p+\\H\ 2 - [H 3 H] = 
, [(« 3 - a)H tg ] - [H 3 u tg ] = 0, [H 3 ] = 0. 
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1. Fast Lax' shocks. Consider an extreme shock. Changing x to — x if 
necessary, the Lax condition read: 

(8.18) I 181 ' 

uj + cj < a. 

In particular, this implies that the front is not characteristic on both side, 
and that the nonglancing conditions in Lemma l8.2l are also satisfied on both 
side and the multiple modes are totally nonglancing. Therefore : 

Proposition 8.5. For extreme Lax shocks, the assumptions of Theorem \l.l\ 

are satisfied. 

2. Slow Lax' shocks. Consider a shock associated to one of the middle 
eigenvalue A± s . Changing x to — x if necessary, the Lax condition read: 

(8.19) % I + , +, 

uj + cj < a < uj + \vj\. 

On both side we have |u3 — o~\ < {v^l, therefore 

Proposition 8.6. For slow Lax shocks, the decoupling condition is never 
satisfied. 



8.5 The H -»• limit 

When H = 0, the system (|8.1|l reduces to isentropic Euler's equations and 
()8.17|) to the corresponding Rankine Hugoniot condition. 
When H = 0, the eigenvalues are 

(8.20) A = A ±1 = A ±2 = u ■ e, A ±3 = A ± c\£\ . 

In particular 8^ 3 Xq = 113. Moreover, at XJ_ = (p,u,0), the tangent character- 
istic polynomial A in (|4.4|) is (r + n-^) 5 . Therefore, if U3 7^ 0, the eigenvalue 
Ao is totally nonglancing. 

Lemma 8.7. Consider a state U_ = (p,u,0). Suppose that 

(8.21) u 3 i {-c,0,+c}. 

Then, for U in a neighborhood ofU_, the boundary x$ = is non characteris- 
tic for the hyperbolic linearized equation and the eigenvalues X±s are simple. 
Moreover, for all £ 7^ ; the multiple eigenvalue Ao is totally nonglancing at 
U. 
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